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1.  INTRODUCTION 


This  Final  Technical  Report  contains  the  summary  of  the  research  results  obtained  under 
the  University  Research  Initiative  (URI)  Grant  AFOSR-90-0712,  entitled  "A  Comprehensive 
Study  of  Matrix  Fracture  Mechanisms  in  Fiber-Reinforced  Ceramic  Matrix  Composites".  The  URI 
Grant  covered  the  period  from  15  February  1990  to  30  April  1993. 

1.1  Research  Objective. 

The  impetus  of  the  research  stems  from  the  Air  Force’s  needs  for  a  class  of  structural 
materials  that  can  function  in  severe  temperature  environments;  and  fiber  reinforced  ceramic  matrix 
composites  have  thus  become  a  prime  subject  for  research  and  development.  Although  some 
ceramics  are  thermally  stable  up  to  a  temperature  beyond  2000  °C,  they  are  invariably  brittle  and 
lack  the  required  strain  capability  as  a  structural  material.  However,  when  it  is  reinforced  with  a 
network  of  equally  thermal  stable,  high  strength  and  (relatively)  high  strain  fibers,  the  brittle  nature 
of  the  ceramic  matrix  may  be  circumvented  and  the  toughness  of  the  composite  increased.  But  the 
toughening  mechanisms  due  to  fiber  reinforcement  can  be  influenced  by  a  host  of  intrinsic  and 
extrinsic  factors,  such  as  composite  processing,  placement  of  the  fibers,  interfacing  with  the 
matrix,  and  the  mechanical  interactions  that  occur  in  the  composite  microstructure  under  the  applied 
loading.  Hence,  there  is  a  need  to  identify  and  characterize  these  influencing  factors;  to  establish 
their  ranges  of  influence;  and  to  include  them  in  a  generic  model  that  can  describe  the  toughening 
mechanisms  and  provide  a  predictive  capability. 

Thus,  the  objective  of  the  present  research  is  to  conduct  a  comprehensive  study  towards 
understanding  the  complex  toughening  mechanisms  in  fiber  reinforced  ceramic  matrix  composites. 
As  the  necessary  first  step,  the  present  study  focuses  on  the  physical  nature  of  matrix  microcracks 
and  the  factors  influencing  their  initiation.  The  scope  of  the  research  thus  entails  composite 
processing,  testing  of  matrix  microcracking,  identification/characterization  of  the  influencing 
factors,  and  the  construction  of  a  predictive  model  suitable  for  numerical  simulation  of  the  matrix 
microcracking  processes  as  observed  in  the  tests.  All  composites  so  studied  are  limited  to 
unidirectionally  continuous  filament-reinforced  ceramic  systems. 

1.2  Research  Approach. 

The  approach  taken  in  the  study  is  one  involving  interactive  efforts  in  experiment, 
theorization  and  simulation.  The  common  thread  linking  these  efforts  is  the  control  and  variation  of 
those  parameters  in  processing,  material  selection  and  loading  which  are  deemed  important  in  the 
matrix  microcracking  processes. 

The  control  and  variation  of  the  influencing  parameters  begins  with  composite  fabrication, 
where  several  material  and  geometric  factors  known  to  influence  matrix  microcracking  are  varied. 
Specifically,  through  the  selection  of  different  fiber-matrix  combinations,  mismatch  of  the  thermal 
and  mechanical  properties  between  fiber  and  matrix  is  controlled;  through  the  selection  of  fibers  of 
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different  sizes  (but  with  the  same  matrix),  local  fiber  geometry  is  altered;  through  the  selection  of 
matrices  of  different  toughness  properties  (but  with  the  same  fiber),  matrix  influence  on  localized 
fracture  is  controlled;  through  treatment  of  the  fiber  surface  prior  to  fabrication,  the  fiber-matrix 
interface  bonding  strength  is  controlled;  and  through  a  unique  processing  method,  the  fiber  volume 
content  in  each  of  the  composite  systems  is  varied  uniformly  from  null  to  about  50%.  Furthermore, 
in  order  to  establish  the  influence  range  of  each  of  the  controlled  parameters,  a  statistical  data  base 
in  each  case  is  generated. 

As  a  result,  the  test  program  becomes  extensive  and  comprehensive.  In  particular,  testing 
and  instrumentation  are  designed  to  emphasize  in-situ  observations  in  real-time;  so  the  actual  matrix 
microcracking  mechanisms  can  be  delineated  visually  and  a  complete  record  for  the  cracking 
process  obtained  as  a  function  of  the  applied  loading.  Generic  model  is  constructed  based  on  the 
observe  microcracking  mechanisms  and  their  relationship  with  respect  to  the  influencing  factors; 
the  successful  model  is  then  applied  to  simulate  each  of  the  test  cases  by  means  of  a  numerical 
procedure  (i.e,  the  finite  element  method).  The  following  provides  a  brief  outline  of  the  various 
tasks  performed  in  the  study. 

1.3  Research  Tasks. 

During  the  three-year  period  of  the  present  research,  numerous  tasks  have  been  performed. 
Key  areas  of  research  conducted  in  the  period  are  briefly  discussed  below: 

In-House  Composite  Fabrication.  This  task  is  the  keystone  of  the  present  research, 
in  that  composite  samples  of  both  quality  and  quantity  must  be  routinely  made  for  testing.  As  is 
mentioned,  selection  of  the  matrix,  fiber  and  interface  treatment  is  so  large  and  variation  of  fiber 
volume  content  in  the  composites  changes  from  null  to  50%,  the  number  of  test  samples  needed  is 
considerable.  During  the  first  two  years  of  the  present  study,  much  of  our  effort  is  devoted  in 
developing  the  fabrication  procedures  to  insure  sample  quality,  in  addition  to  preparing  the  required 
sample  population. 

Testing  and  Instrumentation.  With  the  focus  placed  on  the  "initiation  mechanisms"  of 
matrix  microcracking  in  unidirectional  filament  reinforced  systems,  die  test  condition  chosen  for  all 
samples  is  tension  loading  along  the  fiber  direction.  This  is  achieved  by  3-point  bending  test  of  a 
slender  beam.  The  difficulty  associated  with  the  test  is  that  (a)  to  capture  the  moment  of  microcrack 
initiation;  and  (b)  to  follow  crack  progression  in  the  composite  microstructure  during  loading.  This 
difficulty  is  further  compounded  by  the  fact  that  some  of  the  tests  must  be  conducted  in  high 
temperature  conditions.  As  will  be  discussed  later  in  this  report,  special  test  fixture  and  real-time 
observation  instrumentations  are  developed  during  the  research  period. 

Modeling  and  Simulation.  To  develop  a  theoretical  model  that  can  correctly  link  the 
microstructure,  the  influential  factors  and  the  applied  loading  to  the  physical  events  of  matrix 
microcracking  is  the  central  piece  of  the  present  research.  Since  the  desired  model  must  be  based 
on  mechanics  principles  and  the  actual  cracking  mechanisms,  correlation  between  the  model  and 
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experiment  is  important.  The  correlation  is  provided  by  numerical  simulation  of  the  actual  tests,  so 
that  results  from  the  simulation  can  be  compared  directly  with  experiment;  any  deficiency  of  the 
theoretical  model  is  amended  through  correlative  iterations.  The  value  of  such  a  model  lies  not  only 
in  its  predictive  capability  of  matrix  microcracking,  but  also  in  its  explicit  inclusion  of  the  various 
influencing  factors  which  play  a  role  in  the  toughening  mechanisms  of  the  composites. 

Baseline  Property  Characterization.  In  order  to  provide  the  necessary  input  data  in 
the  numerical  simulation  of  the  tests,  a  considerable  amount  of  effort  is  devoted  to  obtain  the 
relevant  material  properties  for  the  matrices  and  the  fibers  used.  Although  most  of  the  thermoelastic 
constants  for  the  matrices  and  the  fibers  at  room  temperature  can  be  determined  by  routine 
methods,  exceptions  are  abound.  Among  the  most  difficult  ones  are  (a)  the  measurement  of  the 
thermal  expansion  coefficients  of  the  fibers;  and  (b)  the  measurement  of  the  fracture  toughness  of 
the  bulk  matrices  in  high  temperature.  The  former  is  needed  to  assess  the  thermal  residual  stress 
state  of  the  composites,  while  the  later  provides  the  necessary  material  condition  for  matrix 
microcracking. 

Interface  Mechanics  Study.  Fiber-matrix  interfacing  is  a  key  parameter  governing  the 
load-transfer  mechanisms  in  composite  materials.  Proper  mechanics  characterization  of  the 
interface  is  essential  for  the  understanding  of  the  failure  processes  and  the  toughening  mechanisms 
in  composites.  An  effort  is  made  to  rigorously  study  interface  mechanics;  extensive  experiment  and 
theoretical  modeling  are  perform.  Details  on  this  effort  are  presented  later  in  thL  report. 

1.4  The  Research  Team. 

The  following  faculty,  graduate  students  and  research  associates  have  participated  in  the 
research  as  a  team: 

*  A.  S.  D.  Wang  (PI),  Soffa  Professor  of  Mechanical  Engineering 

*  M.W.  Barsoum  (Co-PI),  Associate  Professor  of  Materials  Engineering 

*  T.  M.  Tan  (Co-PI,  departed  Feb.  92),  Associate  Professor  of  Mechanical  Engineering 

*  P.  Kangutkar,  graduate  student  (graduated/Ph.D.  Sept.  91),  Materials  Engineering 

*  X.  G.  Huang,  graduate  student  (graduated/Ph.D.  Feb.  92),  Mechanical  Engineering 

*  P.  V.  Kishore,  graduate  student  (graduated/Ph.D.  Mar.  93),  Mechanical  engineering 

*  T.  Dong,  graduate  student  (Sept.  91 -Mar.  93  with  M.S.),  Materials  Engineering 

*  H.  M.  Chou  (departed,  June.  91),  Post -doctorate  Fellow,  Materials  engineering 

*  J.  C.  Gong  (departed,  Sept.  91),  Post-doctorate  Fellow,  Mechanical  Engineering 

*  A.  Elkind  (Sept.  91 -Mar.  93),  Post-doctorate  Fellow,  Materials  Engineering 

The  team  is  organized  into  five  task-groups,  each  is  assigned  to  specified  responsibility  as 
follows:  (1)  fabrication  and  testing  (Barsoum,  Kangutkar,  Dong);  (2)  thermoelastic  properties  of 
fibers,  matrices  and  composites  (Barsoum,  Kangutkar,  Elkind);  (3)  fabrication  and  testing  of 
interface  properties  (Chou,  Tan,  Gong);  (4)  model  and  simulation  of  interface  mechanics  (Wang, 
Kishore);  and  (5)  modeling,  simulation  and  correlation  of  matrix  microcracking  (Wang,  Huang). 
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A.  S.  D.  Wang  is  the  overall  manager  of  the  research  team. 

1.5  Research  Reporting. 

Throughout  the  course  of  the  present  research,  major  findings  are  reported  via  several 
technical  media:  (a)  AFOSR  Annual  Technical  Reports;  (b)  Ph.D.  Dissertations;  (c)  Archival 
Journal  Papers;  (d)  Proceeding  of  Technical  Symposia;  (e)  Presentation  at  Technical  Conferences; 
and  (f)  Invited  Seminar  at  Universities  and  Research  Institutes.  A  complete  bibliography  resulting 
from  the  present  research  is  contained  in  Section  4. 

In  what  follows,  a  brief  summary  of  the  major  results  obtained  in  the  research  is  given. 
Whenever  possible,  specific  details  are  referred  to  the  appropriate  publications  that  are  listed  in 
Section  4. 


2.  MAJOR  RESULTS  AND  DISCUSSIONS 
2.1.  Composite  Fabrication. 

During  the  first  year  of  this  research,  an  in-house  processing  facility  is  setup  to  fabricate 
unidirectional  continuous  fiber  reinforced  ceramic  matrix  composites  with  the  desired  quality.  The 
processing  procedure  is  then  followed  to  fabricated  numerous  test  samples  for  subsequent  matrix 
microcracking  study.  A  brief  outline  of  the  fabrication  aspects  of  the  research  is  given  below: 

Fabrication  Procedure.  The  composite  fabrication  technique  adopted  in  this  research 
has  been  well  established  in  the  open  literature  (see  Ref.[B-l]  for  specific  details).  For  purpose  of 
discussion,  the  general  procedure  is  as  illustrated  in  Fig.l.  In  brief,  matrix  material  (glass)  in 
powder  form  is  mixed  with  a  polymer  slurry;  the  reinforcing  fiber  is  then  made  to  pass  through  the 
slurry  and  the  slurry  coats  the  fiber  to  form  a  uniform  surface  layer;  the  coated  fiber  is  next 
exposed  to  a  ultraviolet  light  source  which  cures  the  polymer  slurry  coating;  the  coated  fiber  is 
wound  uniformly  onto  a  drum  to  form  a  collimated  prepreg  tape;  from  the  prepreg,  a  green 
unidirectional  laminate  is  made;  the  green  laminate  is  then  heated  in  air  at  500  °C  to  bum  out  the 
polymer  in  the  coating,  leaving  the  glass  powder  adhering  to  the  fiber  surface;  the  laminate  is 
finally  placed  in  a  graphite  die  and  hot-pressed  in  vacuum  to  form  a  consolidated  composite  plate. 

The  fiber-volume  content  in  the  composite  plate  is  determined  by  the  amount  of  glass 
powder  used  in  the  polymer  slurry.  Fiber  contents  from  0  to  50%  are  made  for  each  of  the  selected 
fiber-matrix  combinations.  Composites  of  zero  fiber  content  are  simply  the  bulk  matrix  with  any 
fiber,  hot-pressed  under  the  same  condition  as  those  with  fibers. 

Selection  of  Matrices.  Three  matrix  materials  are  selected  for  fabrication  of  composite 
samples;  and  these  are:  (a)  7740  borosilicate;  (b)  1723  aluminosilicate;  and  (c)  zircon.  All  come  in 
powder  form  and  are  adaptive  to  the  fabrication  procedure  as  outline  before.  The  rationale  for  the 
selection  is  based  on  their  contrasting  material  properties,  especially  the  tensile  modulus,  tensile 
strength,  fracture  toughness  and  the  strain  point  (the  phase  change  temperature  from  liquid  to  solid 
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or  vise  versa);  and  it  is  believed  that  tne  effect  of  these  contrasting  properties  on  matrix  cracking 
should  be  detectable  both  in  experiment  and  in  theory.  The  following  table  lists  some  of  the 
baseline  properties  of  the  three  selected  matrices. 


774Q 

1723 

zircon 

Modulus  (GPa) 

63 

83 

195 

TEC  (10-6/C) 

3.2 

4.6 

4.2 

KIc  (MPaVm) 

0.75 

0.8 

2.5 

Ten.  strength  (MPa) 

80 

35 

265 

Strain  point  (°C) 

520 

650 

900 

Selection  of  Fibers.  Four  fiber  systems  are  selected  for  fabrication  of  composite 
samples  and  they  are  (a)  SCS-6  SiC;  (b)  SCS-9  SiC;  (c)  Nicalon  SiC;  and  (d)  HMU  graphite.  The 
SCS-6  and  SCS-9  fibers  are  monofilaments,  while  the  Nicalon  and  HMU  carbon  fibers  come  in 
tow.  The  following  table  lists  some  of  their  contrasting  properties: 


SCS16 

S£S=2 

Nicalon 

HMU-C 

Modulus  (GPa) 

400 

390 

200 

380 

TEC  (10^/C) 

6.5 

6.5 

3.1 

0.1 

Diameter  (pm) 

75 

35 

8 

4 

Similarly,  the  contrasting  properties  in  these  fibers  are  believed  to  cast  contrasting  effects 
on  matrix  microcracking,  detectable  in  experiment  and  in  theoretical  predictions. 

Interface  Treatment.  To  study  the  effect  of  interface  bonding  on  matrix  microcracking, 
the  SCS-6/7740  borosilicate  composite  is  selected  for  interface  treatment.  The  treatment  procedure 
involves  heating  the  SCS-6  fiber  in  air  at  600°C  for  about  1  minute,  to  achieve  a  thin  oxidized 
surface  layer;  upon  processing  into  composite,  the  oxidized  surface  layer  of  the  fiber  provides  a 
much  stronger  interface  bonding.  By  means  of  the  fiber  pullout  test,  fiber  so  treated  results  in  an 
interface  shear  strength  exceeding  100  MPa,  while  the  untreated  fiber  provides  only  about  8  MPa. 
The  details  of  the  surface  treatment  and  the  subsequent  interface  shear  tests  are  reported  in  Ref. 
[C-2], 

Composite  Systems  Fabricated.  Eight  composite  systems  are  fabricated  during  the 
research  period,  using  the  three  selected  matrix  materials  and  the  four  selected  fibers.  These  are 
listed  as  follows: 

1.  SCS-6/7740  (without  interface  treatment) 

2.  SCS-6/7740  (with  interface  treatment) 
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3.  SCS-9/7740 

4.  Nicalon/7740 

5.  HMU-C/7740 

6.  SCS-6/1723 

7.  Nicalon/1723 

8.  SCS -6/zircon 

In  each  system,  the  fiber  volume  content  is  varied  uniformly  from  0%  to  approximately 
50%.  This  resulted  in  fabricating  some  16  plates  of  4"x4"  in  size  for  each  composite  system  (from 
each  plate,  5  to  6  test  specimens  could  be  made).  More  than  600  test  samples  in  all  are  fabricated 
and  tested  during  the  research  period. 

Detailed  documentation  for  the  processing,  fabrication  and  sample  preparation  is  contained 
in  the  Ph.  D.  thesis  by  Kangutkar  [B-l]. 

2.2.  Matrix  Microcracking  Test  Specimens  and  Procedures. 

Matrix  microcracking  tests  conducted  in  this  research  use  the  3-point-bend  loading  method 
where  matrix  microcracking  events  are  induced  on  the  tensile  face  of  the  3-point  specimen  under 
center  loading.  The  test  specimens  are  in  the  form  of  a  beam  with  the  nominal  dimensions  of  60 
mm  long,  5  mm  wide  and  2  mm  deep.  Specimen  preparation  before  testing  is  a  critical  in  terms  of 
the  quality  of  test  data.  For  instance,  each  specimen  must  be  finely  polished  in  order  to  eliminate 
flaws  and  excess  glass  on  the  surface  and  near  the  comers.  Details  for  the  specimen  preparation  are 
fully  described  in  Ref.  [B-l]. 

Most  of  the  3-point  bend  tests  are  conducted  under  room  temperature;  a  number  of  selected 
composite  systems  are  also  tested  in  higher  temperatures.  In  order  to  monitor  the  matrix  cracking 
events  in  real-time  and  in-situ,  two  different  recording  methods  are  used,  depending  on  the  test 
temperature.  A  brief  outline  of  these  two  methods  is  given  as  follows: 

1.  At  room  temperature,  tests  are  conducted  inside  the  chamber  of  a  scanning  electronic 
microscope  (SEM).  A  specially  designed  test  fixture  is  housed  inside  the  SEM,  with  a  remove 
load-control  devise;  visual  details  of  matrix  microcracking  could  be  observed  in-situ  under  the 
power  of  the  SEM  up  to  lOOOx.  The  test  fixture  used  is  schematically  illustrated  in  Fig.  2. 

2.  At  temperatures  higher  than  room  temperature,  each  test  specimen  is  prepared  with  a  thin 
gold-film  sintered  on  its  tensile  surface  prior  to  testing,  as  illustrated  in  Fig.  3.  Tests  are  then 
conducted  in  a  thermal  chamber  affixed  in  the  loading  train  on  a  MTS  testing  machine.  During  the 
test,  a  devise  is  used  to  measure  the  electric  resistance  of  the  gold  film  in  the  direction  along  the 
beam.  When  matrix  microcracking  occurs  in  the  specimen  and  breaks  the  gold-film,  the  electrical 
resistance  of  the  film  increases  abruptly.  The  corresponding  load  is  recorded  and  used  to  compute 
the  critical  tensile  stress  for  matrix  microcrack  initiation. 

The  SEM  in-situ  test  provides  the  visual  details  of  the  matrix  cracking  mechanisms  at  the 
fiber-matrix-interface  scale.  For  example,  one  of  the  major  findings  in  this  research  is  that  matrix 
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microcracking  initiates  generally  within  a  local  fiber  spacing  as  an  isolated  event;  as  the  applied 
load  increases,  more  such  local  events  occur  within  other  fiber  spacings.  Localized  interface 
debonding  then  forms  near  the  matrix  cracking  sites  at  higher  load  levels;  these  local  matrix 
microcracks  eventually  link  up  through  the  interface  debonds,  resulting  in  a  larger  and  zigzagged 
matrix  crack  across  the  width  of  the  test  sample.  A  schematic  of  the  cracking  initiation  and  linking 
sequence  is  illustrated  in  Fig.  4.  This  finding  is  important  in  determining  the  mechanisms  of  matrix 
cracking  and  in  guiding  the  development  of  the  cracking  simulation  model. 

A  typical  matrix  microcrack  initiation  stress  (MCIS)  data  is  shown  in  Fig.  5  from  testing 
the  SCS-6/7740  specimens  (without  interface  treatment)  in  room  and  higher  temperatures.  The 
MCIS  data  is  plotted  versus  the  fiber  volume  content  Vf  of  the  test  specimens,  with  Vf  varying 
from  0  to  0.5.  Though  the  test  data  is  quite  scattered  they  do  provide  both  a  statistical  basis  and  a 
definitive  behavioral  tend,  useful  for  the  eventual  validauon  of  the  developed  simulation  model. 

Again,  the  experimental  aspects  regarding  the  tests  and  test  data  analysis  are  detailed  in  the 
thesis  of  Kangutkar  [B-l]. 

2.3.  Simulation  Model  for  Matrix  Microcracking. 

The  theoretical  framework  for  modeling  matrix  microcracking  is  first  conceived  in  1988. 
The  original  idea  is  further  developed  during  this  research  due  to  the  extensive  experiment.  Most  of 
the  details  in  the  model  development  are  described  in  the  Ph.  D  thesis  by  Huang,  Ref.[B-2], 

The  developed  model  differs  from  models  of  the  past  in  one  major  aspect.  Namely,  the 
present  model  describes  the  composite  microstructure  at  the  fiber-matrix  interface  scale  and 
accounts  for  the  effect  of  material  defects  both  in  the  matrix  and  on  the  matrix-fiber  interface.  Fig. 
6  shows  the  schematics  of  the  matrix  microcracking  model  and  the  input  par  imeters  included  in  the 
model.  Specifically,  these  parameters  include:  the  thermoelastic  moduli  of  the  matrix  and  the  fiber; 
the  composite  processing  temperature  parameter;  the  fracture  toughness  of  the  matrix;  the 
geometrical  details  of  the  local  microstructure  (e.g.  fiber  diameter  and  fiber  spacing);  and  the 
representation  of  local  defects  both  in  the  matrix  and  on  the  fiber-matrix  interface. 

Clearly,  not  all  of  the  parameters  are  deterministic;  some  are  actually  random  in  nature.  This 
causes  uncertainty  in  the  input  data  into  the  model  and  a  bounding  approach  is  developed  to 
compliment  the  model.  Nevertheless,  the  model  is  still  capable  of  describing  the  complex 
interactions  amongst  the  local  defects  that  are  situated  in  an  evolving  composite  microstructure,  and 
actually  mimicking  the  matrix  crack  process  under  the  applied  loading  as  it  would  in  a  real  test 
situation.  Thus,  the  model  can  provides  a  description  of  the  actual  cracking  process  in  quality  and 
quantity,  although  the  augmentation  of  the  model  must  be  complimented  by  a  bounding  technique 
due  to  the  uncertainty  involved  in  the  input  parameters.  Details  of  these  development  are  found  in 
[B-2]  and  [C-4], 

The  successful  model  is  used  to  perform  two  major  tasks.  The  first  is  a  parametric  study  on 
matrix  microcracking  initiation.  The  study  is  quite  extensive  in  that  it  examines  the  microcracking 
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behavior  as  it  is  influenced  by  all  the  parameters  included  in  the  simulation  model.  The  study  not 
only  establishes  the  value-ranges  for  each  of  the  parameters  but  also  allows  an  examination  of  the 
intricate  interactions  amongst  the  parameters.  The  other  task  concerns  the  simulation  of  the  actual 
test  cases  (about  600  cases)  and  the  prediction  of  the  MCIS.  The  latter  is  used  in  comparison  with 
the  experiment. 

A  typical  theoretical-experimental  comparison  is  shown  in  Fig.  7,  where  the  MCIS  data 
from  testing  and  simulation  of  the  SCS-6/7740,  SCS-6/1723  and  SCS-6/zircon  specimens  are 
shown  plotted  versus  the  fiber  volume  fraction,  Vf.  It  is  seen  that  the  agreement  between  theory 

and  experiment  is  rather  good.  In  fact,  close  comparison  of  MCIS  between  the  predicted  and  the 
tested  results  are  consistently  obtained  from  all  the  8  different  composite  systems,  see  [B-2].  This 
is  rather  remarkable  in  view  of  the  diversity  in  the  microsti  icture  features  exhibited  in  the  differing 
composite  systems. 

2.4  Major  Findings  from  Matrix  Microcracking  Results. 

From  over  600  test  data  and  in  conjunction  with  the  results  from  model  simulations,  a 
number  of  interesting  findings  are  obtained.  A  brief  discussion  on  some  of  the  major  findings  is 
presented  in  the  following: 

Test  Specimen  Surface  PolisMng.  By  finely  polishing  the  test  specimens  before 
testing,  it  helps  to  eliminate  surface  defects  and  excess  matrix;  surface  defect  and  excess  glass  can 
cause  matrix  cracking  prematurely  and  Tesult  in  large  data  scatter.  This  is  illustrated  by  the  data 
shown  in  Fig.  8  where  a  comparison  of  the  MCIS  data  is  made  between  the  polished  and  the 
unpolished  SCS-6/7740  specimens.  It  can  be  readily  seen  that  surface  polishing  reduces  the  scatter 
in  the  MCIS  data  (details  of  specimen  surface  polishing  is  presented  in  [B-l ]).  In  the  discussions 
that  follow,  all  data  presented  are  from  testing  polished  specimens. 

Matrix  Microcracking  Begins  as  Isolated  Cracks  in  Fiber  Spacings.  As  has 
been  discussed  earlier,  the  first  matrix  microcracks  are  found  to  form  inside  the  spacing  of  two 
adjacent  fibers.  In  general,  fiber  spacing  in  a  specimen  is  not  uniform  due  to  processing  variability; 
and  the  first  few  matrix  cracks  occur  almost  always  inside  the  largest  fiber  spacings.  Furthermore, 
the  threshold  stress  (MCIS)  is  found  to  correlate  with  the  value  of  the  fiber  spacing  in  which  the 
crack  occurs.  This  later  is  illustrated  in  Fig.  9  which  displays  the  MCIS  data  from  testing  7 
SCS-6/7740  specimens.  In  the  first  figure  the  MCIS  is  plotted  against  the  average  fiber  spacing  on 
the  tensile  surface  of  the  specimens  (the  averaged  fiber  spacing  is  tantamount  to  »he  averaged  fiber 
volume  content,  Vf),  while  in  the  second  figure,  it  is  plotted  against  the  value  of  the  spacing  >n 
which  the  first  matrix  crack  occurs.  From  the  second  plot,  the  ricnd  that  the  MCIS  decreases  with 
increasing  the  fiber  spacing  is  evident. 

This  finding  indicates  the  importance  of  fiber  spacing  uniformity  during  composite 
processing.  As  an  inherent  property  of  the  system,  fiber  spacing  always  exhibits  a  certain  scatter  in 
its  value  (see,  e.g.  the  7  data  points  in  Fig.  9).  Lack  of  sufficient  uniformity  of  fiber  sp  icing 
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causes  scatter  in  the  MCIS  values.  It  is  possible  that  both  the  MCIS  and  its  scatter  may  be 
increased  by  improving  the  fiber  spacing  uriformity  during  processing.  Further  discussion  uu  this 
subject,  together  with  the  SEM  examination  records,  can  be  found  in  Ref.  [B-2]. 

Small  Fiber  Is  Better  Than  Large  Fiber.  From  both  the  parametric  study  on  fiber 
size  and  the  experiment  with  fiber  size  controlled,  it  becomes  evident  that  systems  with  fibers  of 
smaller  diameter  out  performs  systems  with  fibers  of  larger  diameter  in  terms  of  the  MCIS  if  all  the 
other  parameters  are  kept  the  constant.  Fig.  10  shows  the  MCIS  data  (plotted  against  Vf)  from 
testing  a  group  of  SCS-6/7740  specimens  and  a  group  of  Nicalon/7740  specimens;  the  system  with 
the  Nicalun  fiber  provides  substantially  higher  MCIS  than  the  system  with  the  SCS-6  fiber. 
Thougi  both  fibers  are  basically  made  of  SiC,  the  Nicalon  has  a  diameter  ranging  from  10  to  20 

pm  (with  an  average  of  16  pm)  while  the  SCS-6  has  a  diameter  of  150  pm.  The  theoretical  average 
fiber  spacing,  s  depends  on  the  fiber  diameter,  d  and  the  fiber  volume  fraction  Vf  by  the  relation: 

s=d{  V[(V3jt)/(6Vf)]-l  j.Thus,  for  the  same  Vf,  the  average  fiber  spacing  is  directly  proportional  to 
the  fiber  diameter. 

One  might  also  note  that  the  elastic  moduli  of  the  Nicalon  and  the  SCS-6  are  not  the  same; 
the  Nicalon  is  only  half  as  stiff  in  the  longitudinal  direction  as  that  of  the  SCS-6.  If  the  data  in  Fig. 
10  is  converted  into  composite  strain  instead  of  stress  (see  Fig.  1 1),  the  advantage  exhibited  by  the 
Nicalon  is  even  more  pronounced.  Fig.  12  shows  a  similar  comparison  between  the  SCS-6/1723 
and  the  Nicalon/1723  systems.  Here,  again,  the  Nicalon  out  performs  the  SCS-6  in  terms  of  matrix 
cracking  strain. 

The  comparisons  presented  above,  however,  are  not  strictly  based  on  the  difference  in  fiber 
diameter  alone  since  the  Nicalon  and  the  SCS-6  are  different  fibers  in  other  aspects,  such  as 
elasticity,  interface  bonding  characteristics,  processing  preparation,  etc. 

In  order  to  circumvent  the  difficulty  associated  with  these  differences,  another  experiment 
is  conducted  in  which  systems  made  of  the  SCS-6  and  SCS-9  are  fabricated  and  tested.  The  two 

fibers  are  basically  the  same  in  every  aspects  except  that  the  SCS-9  has  a  diameter  of  70  pm,  while 

the  SCS-6  150  pm.  When  these  fibers  are  coupled  with  the  same  matrix  (the  7740)  and  processed 
under  the  same  condition,  the  effect  of  fiber  size  alone  is  isolated. 

Fig.  13  displays  the  comparison  of  the  MCIS  between  the  SCS-6/7740  and  the  SCS-9/7740 
systems  as  plotted  against  Vf,  both  curves  are  predicted  using  the  developed  simulation  model  (this 
prediction  is  based  on  the  parametric  study  performed  before  the  procurement  of  the  SCS-9  fiber). 
The  prediction  indicates  that  the  SCS-9  fiber  out  performs  the  SCS-6  fiber  in  MCIS. 

Fig  14  displays  the  comparison  of  MCIS  obtained  by  experiment.  In  this  case,  the  MCIS 
data  display  a  considerable  scatter;  the  difference  in  MCIS  between  the  two  systems  can  still  be 
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discerned  in  the  Vf  range  above  30%.  In  this  Vf  range,  the  average  fiber  spacing  in  the  SCS-9 
system  becomes  small  enough  to  effect  an  increased  fiber  constraint  against  matrix  microcracking. 
This  is  not  the  case  for  the  SCS-6  system. 

It  should  also  be  mentioned  that  the  interface  bonding  strength  in  the  SCS-6/7740  system  is 
about  the  same  as  that  in  the  SCS-9/7740  system  (this  will  be  further  discussed  later  in  this  report). 
Thus,  the  comparison  presented  in  Fig.  14  can  be  considered  as  based  on  the  difference  in  fiber  size 
alone.  Overlay  of  Fig.  13  and  Fig.  14  will  show  that  the  model  predicted  MCIS  by  falls  reasonably 
well  inside  the  experimental  data  band. 

Effect  of  Interface  Bonding.  The  reinforcing  fibers  will  provide  a  constraining  effect 
on  matrix  microcracking  only  if  the  fiber-matrix  interface  is  adequately  bonded.  If  the  interface 
bonding  is  weak,  the  threshold  stress  for  matrix  microcracking  (MCIS)  is  generally  lowered.  A 
physical  explanation  with  mechanics  quantification  of  the  effect  is  fully  described  in  the  simulation 
model,  see  Ref.  [B-2].  An  experiment  is  also  conducted  to  validate  the  effect  of  interface  on  the 
MCIS.  In  this  case,  composites  made  of  surfaced  treated  and  untreated  SCS-6  fibers  with  the  7740 
matrix  are  fabricated  and  tested  (details  for  fiber  surface  treatment  are  presented  in  [B-l  ];  also,  in  a 
separate  test  involving  single  fiber  pullout,  the  interface  shear  strength  of  the  untreated  fiber  is 

found:  x  =  6  -  8  MPa,  while  that  of  the  surface  treated  fiber:  x  >  100  MPa).  Fig.  15  shows  the 
comparison  of  the  measured  MCIS  from  the  group  with  untreated  fibers  and  the  group  with  surface 
treated  fibers;  the  group  with  surface  treated  fibers  has  substantially  higher  MCIS  than  the  group 
without  surface  treatment. 

Effect  of  Matrix  Flaw  Size.  The  fundamental  notion  that  microflaws  are  present  in  the 
matrix  phase  of  the  composite  and  they  are  the  source  of  matrix  microcracking  is  the  basis  on 
which  the  crack  initiation  mechanism  is  modeled.  Within  that  context,  the  actual  identity  of  the 
microflaws  (size,  shape  and  distribution)  cannot  be  realistically  determined.  However,  as  an 
effective  quantity,  the  microflaws  may  be  considered  as  being  inherent  to  the  matrix  material,  the 
interactions  with  the  fibers  during  processing,  post-processing  handling,  etc.  In  that  sense,  they 
can  come  as  a  random  variable.  It  is  then  argued  that  the  effective  size  of  the  flaws  must  be 
bounded  either  by  (a)  the  size  of  fiber  spacing,  s  and  (b)  the  size  of  the  largest  effective  flaw  in  the 
bulk  of  the  matrix,  a, whichever  is  smaller.  The  former  can  be  evaluated  by  knowing  the  fiber  size 
and  fiber  volume  content;  while  the  latter  can  be  inferred  by  applying  the  classical  fracture 

mechanics  for  brittle  crack,  assuming  that  the  matrix  tensile  strength  om  is  determined  by  the 
matrix  fracture  (mode-I)  toughness  Kjc  and  the  effective  flaw  size  a. 

For  instance,  between  the  7740  and  the  1723  glasses  in  bulk,  the  corresponding  Kjc  are 

0.75  MPaVm  and  0.8  MPaVm,  respectively;  the  corresponding  om  are  80  MPa  and  35  MPa, 


respectively.  Thus,  the  corresponding  dominant  effective  flaw  sizes  are  60  pm  and  300  pm, 
respectively.  This  is  why  the  1723  glass  is  weaker  than  the  7740  glass  in  terms  of  tensile  crack 
resistance.  Hence,  if  reinforced  with  the  same  amount  of  fiber,  the  threshold  stresses  at  matrix 
microcracking  initiation  can  be  elevated  more  significantly  in  the  1723  system  than  in  the  7740 
system. 

Fig.  16  shows  the  comparison  of  the  MCIS  as  a  function  of  fiber  volume  content  between 
the  SCS-6/7740  and  the  SCS-6/1723  specimens.  The  MCIS  is  increased  much  more  dramatically 
in  the  1723  system  due  to  increase  of  fiber  reinforcement  than  in  the  7740  system.  According  to  the 
effective  flaw  concept,  the  average  fiber  spacing  in  the  7740  system  is  larger  than  the  maximum 
effective  flaw  size  in  the  bulk  matrix  when  the  fiber  volume  content  Vf  is  less  than  0.3.  In  that 
case,  the  reinforcing  fibers  are  ineffective  in  constraining  the  initiation  of  matrix  microcracking.  A 
full  discussion  on  this  subject  is  presented  in  Ref.[D-3]. 

2.5  Fiber  Pullout  Tests  and  Interface  Study. 

One  of  the  main  tasks  performed  in  this  research  is  to  provide  the  necessary  baseline 
properties  in  order  for  model  simulation.  These  properties  include  the  thermoelastic  constants  of 
the  matrix  and  fiber  material  used,  the  tensile  strength  properties  of  the  matrix  including  the  mode -I 
fracture  toughness,  and  the  interface  properties  between  the  matrix  and  fiber.  Most  of  these 
properties  can  be  gathered  either  from  prior  studies  or  from  routine  characterization  tests.  The 
interface  property,  however,  is  unique  to  the  particular  composite  processing  in  question  and  the 
interface  property  characterization  procedure  is  by  no  means  routine.  In  fact,  a  significant  effort  is 
made  during  this  research  towards  the  interface  study;  both  experiment  and  model  simulation  are 
conducted.  Some  key  point  in  this  effort  will  be  briefly  discussed  below. 

The  Single  Fiber  Pullout  Test.  The  interface  study  test  program  is  carried  out  using 
the  SCS-6/7740,  SCS-9/7740  and  the  Nicalon/7740  systems.  It  involves  fabrication  and  testing  of 
a  large  number  of  single-fiber  pullout  specimens.  The  objective  of  the  test  program  is  to  examine 
the  underlying  mechanisms  in  the  fiber  pullout  process  as  well  as  to  determine  the  actual  values  of 

the  interface  shear  strength,  x  in  all  the  composite  systems.  The  former  provides  the  physical  basis 
for  the  development  of  a  fiber  pullout  simulation  model,  while  the  latter  is  used  as  input  into  the 
matrix  microcracking  simulation  model. 

Fig.  17  shows  the  schematics  of  the  pullout  test  and  Fig.  18  is  a  photograph  of  the  actual 
test  setup.  Test  results  of  the  interface  shear  strength  and  related  discussions  are  found  in  Refs. 
[A-2J,  [C-l],  [C-2]  and  [H-1J. 

Three  different  mechanisms  which  govern  the  fiber  pulling  out  motion  are  identified:  (a) 
purely  elastic  response  when  the  interface  is  intact  or  undamaged;  (b)  interface  debonding  near  the 
matrix  surface  and  propagating  down  the  fiber's  embedded  length;  and  (c)  frictional  sliding  when 
the  interface  is  debonded.  The  mechanisms  (b)  and  (c)  may  occur  separately  or  concurrently 


depending  on  a  number  of  factors,  such  as  the  residual  stresses  due  to  processing,  the  ambient  test 
temperature,  etc. 

Simulation  of  Fiber  Pullout.  A  simulation  model  is  developed  based  on  the  three 
mechanisms  identified  in  the  fiber  pullout  tests.  The  model  is  then  used  to  sort  out  the  effects  of  the 
various  influencing  factors  governing  the  property  of  the  interface. 

The  simulation  model  is  formulated  on  basis  of  three-dimensional  elasticity;  and  the 
resulting  governing  equations  are  solved  using  a  'solution-matching'  procedure.  This  procedure  is 
new  to  the  computational  mechanics  literature  and  it  offers  both  accuracy  and  efficiency  in 
simulating  the  fiber  pullout  motion.  Briefly,  the  procedure  involves  the  determination  of  the 
stresses  of  the  singular  field  neat  the  fiber  and  the  stresses  in  the  far-field  where  loading  and 
boundary  conditions  are  specified;  accuracy  of  the  computed  stresses  and  the  corresponding 
fracture  quantities  are  assured  by  matching  the  two  solutions  perfectly  within  the  domain  where 
they  overlap. 

The  developed  model  and  the  computational  scheme  are  coded  on  a  high-speed  computer  to 
mimic  the  pullout  test  with  the  mechanisms  that  are  identified  in  the  tests,  including  fiber  pulling 
with  perfect  interface  bonding,  interface  debonding  initiation  and  propagation  and  frictional  sliding 
after  fiber  debonding. 

The  above  development  forms  the  basis  of  the  Ph.  D.  thesis  by  Kishore  Pochiraju  [B-3];  a 
copy  of  which  has  been  submitted  to  AFOSR  as  a  part  of  our  routine  reporting.  A  series  of  archival 
papers  resulting  from  this  research  have  either  appeared  or  been  submitted  for  publication  in  the 
open  literature.  A  copy  of  the  paper  [D-4]  which  details  the  stress-fields  matching  technique  is 
included  in  the  Appendix. 

2.6  Thermal  Expansion  Coefficient  of  the  SCS-6  Fiber.  Thermal  residual  stresses  in 
composites  play  a  significant  role  in  the  formation  of  microscopic  damages  in  the  composite.  In 
order  to  quantify  the  contribution  of  the  thermal  effect  to  matrix  microcracking,  determination  of 
the  thermal  expansion  coefficients  of  the  matrix  and  the  fiber  is  the  essential  first  step.  Generally, 
thermal  coefficients  of  the  commonly  used  materials  can  be  found  either  from  the  open  literature,  or 
from  the  material  supplier.  However,  the  reliability  of  such  sources  cannot  be  certain.  This  has 
been  the  case  in  the  present  research  when  the  SCS-6  fiber  is  used. 

Specifically,  the  longitudinal  thermal  expansion  coefficient  of  the  SCS-6  filament  has  a 
nominal  value  in  the  range  of  4  to  5.5xlO*fy°C,  according  reports  in  the  literature.  The  wide  range 
of  the  reported  values  reflects  the  uncertainty  of  the  source;  and  it  causes  difficulty  in  their  use  in 
order  to  insure  the  computed  thermal  residual  stresses  in  the  composite  with  the  required  accuracy. 
Thus,  an  effort  is  made  in  this  research  to  independently  measure  the  thermal  expansion  coefficient 
of  the  SCS-6  fiber  in  the  range  of  25  to  900  °C.  The  experiment  itself  and  subsequent  follow-up 
investigations  have  uncovered  a  number  of  important  physical  and  chemical  characteristics  of  the 
SCS-6  fiber  in  high  temperature  environments.  In  particular,  the  thermal  expansion  coefficient 
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measured  in  this  research  has  a  value  of  6.5±  0.5  xlO'fyoc,  much  higher  than  any  of  the  reported 
values.  Through  a  number  of  direct  and  indirect  checks,  it  is  concluded  that  the  higher  value 
measured  is  probably  correct;  the  reason  for  the  high  value  can  be  attributed  to  the  formation  of 

P-SiC  grain  growth  during  high  temperature  heating,  which  subsequently  changes  the  overall 
physical  property  of  the  fiber. 

The  details  of  the  experiment  and  related  discussions  on  the  results  are  contained  in  Refs. 
[C-5],  [D-5]  and  [D-6].  A  copy  of  Ref.[C-5]  is  included  in  the  Appendix. 

3.  CONCLUSIONS 

In  the  three  years  of  the  present  research,  we  have  focussed  on  achieving  the  main  objective 
set  out  at  the  beginning  of  the  program,  which  is  to  first  understand  and  then  describe  the  complex 
fracture  toughening  mechanisms  in  fiber-reinforced  ceramic  matrix  composites.  In  this  context, 
matrix  microcracking  is  singled  out  for  study  because  it  is  the  fundamental  form  of  damage  in  such 
composites.  It  is  necessary  to  plan  and  to  execute  an  experimental  program  to  identify  the  key 
factors  affecting  matrix  microcracking  and  to  develop  a  simulation  model  that  describes  the  physics 
of  the  cracking  phenomenon.  In  the  experiment,  control  of  the  influencing  factors  in  the  test 
specimens  is  essential  to  assess  their  effects  on  matrix  microcracking;  it  leads  to  an  extensive 
specimen  fabrication  plan.  At  the  same  time,  a  visual  (physical)  comprehension  of  the  cracking 
mechanisms  at  the  fiber-matrix  interface  level  is  helpful;  and  it  leads  to  painstaking  experimental 
procedures  in  order  to  observe  and  follow  the  cracking  process  in- situ  and  in  real-time. 

As  for  the  development  of  the  simulation  model,  the  representation  of  the  composite 
microstructure  and  the  inclusion  of  the  observed  cracking  mechanisms  and  the  factors  controlling 
matrix  microcracking  result  in  a  complex  micromechanics  model  which  requires  (a)  extensive  input 
data  which  must  be  complete  and  accurate;  and  (b)  a  suitable  numerical  routine  that  can  be  executed 
on  the  computing  machine. 

It  is  important  to  point  out  that  some  of  the  necessary  input  data  are  in-situ  properties  and 
are  not  easily  determined  (such  as  the  interface  property).  In  addition,  the  basic  assumptions  based 
on  which  the  simulation  model  is  developed  cannot  be  all  validated  through  experiment.  Thus,  the 
adequacy  of  the  matrix  microcracking  model  can  be  judged  only  based  on  test  data  collected  in  the 
program.  In  that  sense,  it  may  be  said  that  comparison  between  the  test  data  and  the  model 
simulation  is  quite  good  and  extensive. 

Two  important  areas  remain  yet  to  be  explored  more  fully.  One  area  is  the  study  of  the 
linking  mechanisms  of  the  isolated  matrix  cracks;  crack-linking  is  a  major  energy  dissipation 
mechanism.  The  other  area  is  the  study  of  high  temperature  on  the  cracking  mechanisms.  Both 
areas  are  touched  upon  lightly  in  this  research;  they  remain  to  be  followed  up  in  the  future. 
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Fig.  1 .  Schematic  flaw-chart  for  composite  fabrication 
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Fig.  2  Schematics  of  the  SEM  3-point  bend  test  fixture  with  specimen 
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Fig.  3.  Schematics  of  the  test  specimen  with  a  gold-film  for  surace  crack  detection. 
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Fig.  4.  Schematics  of  the  matrix  microcracking  sequence  as  seen  in  the  SEM. 
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Fig.  5.  A  collection  and  presentation  of  the  matrix  cracking  initiation  data  (SCS-6/7740). 
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Parameters: 

constituents  elastic  moduli  -  (E,  v)f  m 
strength  properties  -  ym 

thermal/processing  -  Aa,  AT 
microstructure  -  a  probabilistic  local  unit-cell: 

-  Vf,  R  (fiber  radius),  s  (fiber  spacing) 

-  microflaws  in  matrix  and  interface  (a,  b) 


Fig.  6.  The  unit-cell  model  for  matrix  microcrack  simulation 


Fig.  7.  Comparison  between  the  predicted  and  experiment  MCIS  results 
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Fig.  9.  Effect  of  fiber  spacing  on  MCIS.  (above)  MCIS  plotted  vs.  average  spacing; 
(below)  MQS  plotted  vs.  spacing  in  which  crack  oc erred 
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Fig.  12.  Effect  of  fiber  diameter  on  matrix  cracking  strain.  SCS-6  vs.  Nicalon  with  1723 
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Matrix  cracking  stress,  MPa 


Fig.  14.  Effect  of  fiber  diameter  on  MCIS  as  illustrated  by  test  results. 
SCS-6  and  SCS-9  with  7740 
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Fig.  15.  Effect  of  Fiber  surface  treatment  on  MSIC. 
SCS-6  with  7740 


Fig.  18.  Machine  set-up  for  the  single  fiber  pullout  test 
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APPENDIX 


The  following  papers  are  appended  to  this  report  for  purpose  of  reference  (ail  papers  listed  above 
may  be  furnished  upon  request): 


1 .  "A  local-global  matching  method  for  the  single  fiber  pullout  problem  with  perfectly  bonded 
interface.''  by  K.  Pochiraju,  A.  C.  W.  Lau  and  A.  S.  D.  Wang,  J.  Computational  Mechanics 
(accepted). 

2.  "Thermal  expansion  of  silicon  carbide  monofilaments  and  silicon  carbide  borosilicate 
composites."  bv  A.  Elkind.  M.  Barsoum  and  P.  Kangutkar.  J.  Am.  Ceramic  Soc.  Vol.  75, 
1992.  pp.287 1-73. 
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A  Local-Global  Matching  Method  for  the  Single  Fiber 
Pullout  Problem  with  Perfectly  Bonded  Interface 
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'•’Center  for  Composite  Materials,  201,  CMSL 
University  of  Delaware, 

Newark,  DE  19716  USA 

*  Corresponding  Author 

§Department  of  Mechanical  Engineering  and  Mechanics 
Drexel  University, 

Philadelphia,  PA  19104  USA 

Abstract 

This  paper  presents  a  local-global  matching  method  to  effectively  solve  both  the 
detail  spatial  structure  of  the  locally  singular  stress  field  and  the  complete  global  stress 
distribution  in  the  single  fiber  pullout  model.  The  main  motivation  to  solve  for  the 
locally  accentuated  stress  field  is  the  belief  that  locally  elevated  stresses,  strains  and 
energy  are  likely  to  induce  damage  formation  and  growth.  The  method  consists  of  three 
components:  a  local  analysis,  a  global  analysis,  and  proper  local-global  matching  of  the 
local  asymptotic  field  to  the  global  complete  stress  field.  The  method  is  applicable  for 
various  fiber-matrix  interfacial  conditions:  perfectly  bonded,  partially  debonded  with 
interfacial  crack,  or  debonded  with  frictional  interfacial  sliding.  Results  presented  here 
are  those  of  perfectly  bonded  fiber-matrix  interface.  For  tiiis  problem,  the  local  stress 
field  is  asymptotically  singular  (unbounded)  at  the  immediate  vicinity  where  the  fiber 
protrudes  from  the  matrix.  In  the  local  analysis,  use  of  complex  potential  and  asymptotic 
expansion  shows  that  when  the  fiber  is  stiffer  than  the  matrix  the  most  dominant  singular 
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stresses  at  the  fiber  protrusion  point  are  of  the  form:  atJ{p,9)  =  K  pxdij(9)  where  the 
exponent  X  is  real-valued  and  X  <  0.  The  local  analysis  determines  the  analytic  spatial 
structure  of  the  local  field:  the  radial  dependence  px  and  angular  variations  of  the 
stresses  <^(0).  The  actual  magnitude  of  the  local  stress  field  is  scaled  by  the  amplitude 

factor  K  which  depends  upon  externally  applied  load  and  overall  boundary  conditions. 
The  global  analysis,  performed  by  a  finite  element  model,  can  be  subjected  to  arbitrary 
fiber-pulling  load  and/or  thermal  load.  The  local-global  matching  method  is  based  on 
matching  the  angular  variation  of  the  local  asymptotic  field  to  the  finite-element 
computed  full-field  stresses  to  determine  K.  Several  monitoring  parameters  are 

developed  to  quantify  the  quality  of  the  matching.  In  many  composite  systems,  there  are 
two  singular  terms  in  the  local  field:  a  =  Kpx'&  (A;0)  +  Kpx'd2(X:9). 

Consequently,  local-global  matching  procedures  for  one-term  matching  and  two-term 
matching  are  developed  in  this  paper.  Then,  the  method  is  generalized  to  determine 
complex-valued  K  for  stress  field  having  complex-valued  X.  It  is  pointed  out  that  the 
local-global  matching  method  also  applies  to  problems  with  material  nonlinearity. 
Numerical  examples  are  presented  to  illustrate  germane  features  of  the  method. 

1  Introduction 

The  load  transfer  between  fiber  and  matrix  is  of  importance  to  the  evolution  of 
damage  in  composites.  In  particular,  the  fiber-matrix  load  transfer  during  fiber  pullout, 
the  interfacial  damage  mechanisms,  the  interface  strength,  and  the  overall  toughness  of 
brittle  matrix  composites  are  intimately  related  (Kerans,  Hay,  Pagano  and  Parthasarathy 
(1989)  and  Evans  and  Marshall  (1989)).  The  pulling-out  (or  pushing-in)  of  a  single  fiber 
embedded  in  a  concentric  cylinder  of  matrix  has  been  used  to  study  the  fiber-matrix  load 
transfer  and  to  characterize  the  fiber-matrix  interfacial  strength.  For  the  single  fiber 
pullout  model,  this  study  presents  an  effective  method  to  solve  for  not  just  the  overall 
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stress  distribution  but  also  the  detail  analytical  form  of  accentuated  local  stresses  at 
geometric  and  material  discontinuity  where  damage  is  most  likely  to  develop. 

The  problem  of  a  single  fiber  embedded  in  a  matrix  was  examined  by  Greszczuk 
(1969).  McCartney  (1989),  Hutchinson  and  Jenson  (1990)  and  Kurtz  and  Pagano  (1991), 
among  others.  In  these  solutions,  the  emphasis  was  not  on  the  detail  stress  distribution  at 
the  local  regions  where  damage  may  initiate  or  propagate,  hence  the  exact  nature  of  the 
local  accentuated  stress  field  near  the  fiber  protrusion  point  (one  of  the  local  regions)  was 
not  rigorously  determined.  In  our  work,  closed-form  expressions  were  determined  to 
solve  for  the  radial  and  angular  variations  of  locally  singular  stresses  for  the  single  fiber 
pullout  problem.  Unlike  crack-tip  solutions,  in  this  problem  two  singular  stress  terms 
with  different  radial  dependence  are  found  for  many  fiber  and  matrix  combinations.  This 
paper  focuses  on  developments  relevant  to  the  local-global  matching  method.  Other 
details  on  the  analyses  and  characteristics  of  the  stress  fields  are  presented  in  Pochiraju 
(1993)  and  Pochiraju,  Lau  and  Wang  (1992,  1993). 

The  general  approach  of  solution  involves  three  essential  steps.  The  first  step 
concerns  the  detail  form  of  local  stress  field  in  the  singular  region.  In  the  local  analysis,  a 
power-type  asymptotic  eigenvalue  expansion  of  the  Muskhelishvili-Kolosov  stress 
potential  functions  is  utilized  to  solve  for  the  singular  local  stress  field.  The  local 
analysis  determines  the  exponent  of  the  singularity  and  the  associated  spatial  variation  of 
the  dominant  stress  field.  The  second  step  is  a  global  analysis  using  finite  element 
method  to  determine  the  complete  stress  distribution  over  the  entire  domain 
corresponding  to  any  specific  mechanical  and/or  thermal  loading  and  boundary 
constraints.  In  the  third  step,  a  local-global  matching  method  based  on  the  angular 
variation  of  stresses  is  performed.  A  zone  is  determined  near  the  singular  point  in  which 
the  angular  variations  of  the  finite  element  computed  stresses  match  perfectly  with  those 
of  the  local  asymptotic  singular  field.  The  perfect  match  determines  the  amplitude  factor 
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K  of  the  local  singular  stress  field,  thus  giving  completely  the  analytic  local  stress 
distribution. 

This  local-global  matching  method  is  versatile  and  can  be  used  to  solve  the  local 
singular  stress  fields  in  a  variety  of  problems.  This  analytical-computational  method  is 
easily  applicable  for  arbitary  mechanical  and/or  thermal  loading,  and  for  a  variety  of 
fiber-matrix  interfacial  conditions.  Results  of  applying  the  local-global  matching  method 
to  frictional  sliding  and  partially  cracked  fiber- matrix  interface  problems  are  reported  in 
Pochiraju  (1993).  With  minor  modification,  the  local-global  matching  method  presented 
here  is  also  applicable  to  material  nonlinear  problems.  One  of  the  authors  (Lau)  has 
applied  local-global  matching  to  power-law  creeping  alloys  with  sliding  grains  (Lau, 
Argon  and  McClintock  (1983))  and  power-law  plastically  deforming  metal  matrix 
composites  (Tan,  Lau  and  Rahman  (1989)).  The  theoretical  basis  and  recent 
developments  of  the  local-global  matching  method  are  presented  in  Sections  3,  4  and  5. 
As  a  starting  point.  Section  2  summarizes  the  determination  of  local  singular  stress  field 
in  the  fiber  pullout  problem. 


2  Local  Analysis 

The  physical  model  considered  is  shown  in  Figure  1.  A  single  fiber  of  radius  R0 
is  embedded  in  a  concentric  cylinder  of  matrix  material  of  outer  radius  -^==,  where  Vf  is 

the  fiber  volume  fraction  of  the  composite  material.  A  global  cylindrical  coordinate 
system  {  R,  Z,  0},  with  Z  along  the  axis  of  the  fiber  is  adopted  to  describe  the  overall 
geometry.  The  view  at  right  shows  a  cross-section  of  the  meridian  plane  A-A  at  any 
specific  <{>  value.  As  shown  in  the  cross-sectional  view  A-A,  the  embedded  length  of  the 
fiber  is  denoted  by  t,  and  its  protruded  length  is  assumed  to  be  10Ro  or  longer.  To 

eliminate  rigid  body  motion,  let  the  Z  =  t  plane  of  the  matrix  cylinder  be  constrained  not 
to  move  in  the  Z-direction.  If  the  protruded  fiber  is  pulled  outward  (or  pushed  inward), 
and/or  the  composite  is  subjected  to  a  thermal  load  AT.  load  transfer  from  the  fiber  to  the 
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matrix  takes  place  through  their  mutual  interface  at  R  =  R0.  Since  the  geometry  and 
loading  in  this  model  has  global  axisymmetry,  we  only  need  to  solve  for  the  solution  in 
the  meridian  plane  of  any  arbitrary  <$>. 

The  point  where  the  fiber  protrudes  from  (or  enters  into)  the  matrix  is  labeled  as 
O.  The  fiber  protrusion  point  O  is  a  site  of  geometric  and  material  discontinuity  where 
local  accentuation  of  stresses  and  strains  occur.  In  the  local  analysis,  our  goal  is  to 
determine  the  detail  spatial  structure  of  the  stress  field  at  the  immediate  vicinity  of  O.  The 
motivation  to  solve  for  the  locally  accentuated  stress/strain  field  is  the  belief  that  locally 
elevated  stresses,  strains  and  energy  are  likely  to  induce  damage  formation  and  growth. 
To  solve  for  this  locally  accentuated  (in  fact,  singular)  stress  field,  it  is  convenient  to 
define  a  local  coordinate  system  { p,  0 }  with  origin  at  point  O. 

In  axisymmetric  problems,  at  locations  where  the  length  scale  is  much  smaller 
than  the  local  curvature,  the  dominant  terms  of  the  governing  equation  in  the 
axisymmetric  model  are  identical  to  the  plane  strain  governing  equation  (  see  equations 
(11)  and  (12)  of  Zak  (1964)  and  section  3.4  of  Atkinson,  Avila,  Betz  and  Smelser 

(1982)).  Within  a  region  —  «  1  around  the  fiber  protrusion  point  O,  the  axisymmetric 

R0 

local  stress  field  can  be  solved  by  treating  the  material  as  locally  under  plane  strain 
condition,  and  use  analytic  methods  originally  developed  for  plane  strain  governing 
equations. 

In  the  local  analysis,  we  applied  Muskhelishvili-  Kolosov  complex  potential 
theory  to  both  the  fiber  and  matrix  domains  (see  Section  A-A,  Figure  1).  In  each  domain, 
the  components  of  stresses  are  expressed  in  terms  of  the  spatial  derivatives  of  two 
complex  potential  functions,  £2(z)  and  (0{z)  (England  (1971),  pp.  53-55).  The  fiber- 
matrix  interfacial  conditions  connect  the  solutions  of  the  two  domains  together. 

A  power-type  solution  is  sought  using  Taylor  series  expansion  of  the  complex 
potential  functions  at  the  immediate  vicinity  of  the  fiber  protrusion  point  O: 

5 


x 


Pochiraju,  Lau  and  Wang 


°(*>  =  1(^+4***) 

_  ■  o) 

6)(z)  =  ^(fin2^  +fl>Sn) 

n=l 

The  normal  and  shear  stress  components  and  the  displacement  components  with  a 
single  term  (n=l,  and  Sj=S  )  expansion  in  (1)  are: 


Te(p,  0 )  =  <J6e(p,  0 )  -  iope(p,  0) 

=  Sps'l^Ae‘iS'u0  +  A*Se~,<s'"0  +  B*e'i,s+1)0)  +  $>5"1(aY<S"i,®+  ASe~‘^~u0  +  Be~‘s+U0  ) 


=  sps-,d1(0)+  sp'-le2(0) 

GU(P .  0)  =  G(ue(p ,  6)  +  iup(p,  0)) 

=  ps(AKe,,s-,l0-^Se-‘‘5-,,0-Fe->,s*u0)  +ps(Ar  e‘rs',>0-  ASe-rs-n0  -  Be-'rs+l)S ) 

Denoting  X  =S-1,  for  a  single  term  in  the  expansion  of  (1): 

ajp,0)ocpka,(d) 


(2a) 


(2b) 

(3) 


Equation  (3)  establishes  the  spatial  characteristics  of  the  stress  fields  in  the  local  regions. 
The  stresses  in  the  local  region  have  a  power-type  dependence  on  the  radial  coordinate  p. 
When  X  <  0,  the  local  stresses  have  a  power-type  singularity.  In  our  analysis,  both  real- 
valued  and  complex-value  stress  singularities  are  allowed  ( X  can  be  real  or  complex). 


The  complex-valued  constants  A,  B,  A  ,  B  ,  their  conjugates,  and  possibly 

complex-valued  X  are  determined  using  the  local  boundary  conditions  and  fiber-matrix 

interfacial  conditions  (see  Section  A- A,  Figure  1).  For  a  perfecdy  bonded  fiber-matrix 
interface,  the  stress  components  <7qq  and  OpQ  ,  and  the  displacement  components  Uq 

and  Up  ,  are  continuous  across  the  fiber-matrix  interface.  The  other  surfaces  of  the 

matrix  and  the  protruding  fiber  are  traction  free.  Denoting  quantities  in  the  fiber  domain 
with  a  superscript  / ,  and  quantities  in  the  matrix  domain  with  superscript  m,  the  local 
interfacial  and  boundary  conditions  can  be  expressed  as: 
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Fiber  -  matrix  interface  conditions: 

At  0  =  0:  Tf#(0  =  O)-r^(0  =  O)  =  O 

Ul(6  =  0)-U™(d  =  0)  =  0 

(4) 

Traction  free  boundary  conditions:  v 

At  6  =  -—:  tT(0  =  — -— )  —  0 

At  6  =  n:  tfe(6  =  n)  =  0 

Substituting  (2)  in  (4),  four  homogenous  equations  with  complex-valued 
coefficients  (eight  real-valued  equations)  can  be  obtained.  These  four  equations  along 
with  their  four  complex  conjugates  form  the  system  of  governing  equations: 

[C]{X}  =  {0}  (5) 

where  {*7=^,  Tf  ,  Bf  ,  B)  ,  Am  ,  A’m,  Bm,  B-;}  (6) 

In  this  eigenvalue  problem,  nontrivial  solutions  exists  when  the  determinant  of  [C] 
vanishes.  The  determinant  of  [C]  has  the  analytically  form  (Pochiraju  (1993)): 

D(cr,/J;S)  =  (P2  -  l)cos(3S/r)  +  {4S 2{P  -  a){P  + 1)  -  2 (p2  -  cr)}cos(2S;r) 

-  (a2  +  p2  +  2a)cos(S;r)  -  2S2{(a  -  P)2  +  (0  + 1)2 }  +  2 p2  +  a2  + 1  (7) 

=  0 

In  the  above  equation  a  and  are  the  Dundurs  parameters  (Dundurs  (1969)) : 

t(k{  -f  l)  —  ( Km  + 1)  _  r(xf  -l)-(vm  -1) 

'K  +  O  +  K+O  ’  r(K;f+l)+(vm  +  l) 

Q 

in  which  r  =  — ,  Kf=3-4vf,  Km  =  3  -  4  vm  ,  and  Gf  are  the  shear  moduli  of 

the  matrix  and  fiber  respectively,  and  vm  and  Vf  are  the  Poisson's  ratios  of  the  matrix  and 
fiber  respectively. 
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Our  interest  is  on  the  dominant  terms  of  (7)  that  give  rise  to  locally  accentuated 
stresses.  Among  the  many  eigenvalues  Sn  of  (7),  the  ones  in  the  range  of  0  <  Re(Sn)  < 

1  produce  singular  stress  field  but  finite  strain  energy  at  the  fiber  protrusion  point  O. 
This  corresponds  to  an  admissible  range  of  -1  <  Re(A)  <  0  for  k  (X  =  S-l).  The  zeros  of 
the  determinant  (7)  are  searched  numerically  using  modified  Powell  algorithm  (Mor£, 
1980)  with  real  and  imaginary  components  of  k  as  independent  variables,  within  the 
admissible  range  of  -1  <  Re(X)  <  0. 

After  the  determination  of  the  eigenvalue  k,  the  eignevector  X  can  be 

determined.  Once  the  complex-valued  coefficients 
Af ,  A'f  ,  Bf  ,  B'f  ,  Am  ,  A'm  ,  Bm  ,  B‘m  in  X  are  determined,  the  angular 

variations  (7,(0)  and  d2(d)  can  be  computed  using  (2).  When  the  fiber  is  stiffer  than  the 
matrix,  k  is  real-valued.  In  this  case,  the  rank  of  the  system  of  governing  equations  (5)  is 
such  that  these  coefficients  can  be  determined  to  within  one  arbitrary  real-valued 
constant.  It  follows  that  the  stresses  are  determinable  to  within  one  arbitrary  amplitude 
factor  K,  which  scales  the  magnitude  of  the  stresses.  The  singular  stress  field 
corresponding  to  Xq  is 

<j'°cal(p,e)  =  K,pK'd\ax\0)  (9) 

To  uniquely  define  this  arbitrary  scaling  factor  K,  we  set  the  angular  variation  of  the 
opening  stress  at  the  fiber-matrix  interface  to  be  unity  ,  that  is  <rM(0  =  0)  =1.  The 
angular  variation  of  the  stresses,  d1  (A,;0),  are  determined  by  the  eigenvalue  A.)  and  the 

coefficients  of  the  eigenvector  X  .  The  magnitude  of  dl  (A,;0)  are  bounded  functions 
of  0.  Solution  of  a  typical  a1  (A,;0)  from  local  asymptotic  analysis  is  depicted  by  the 

lines  in  Figure  4. 

The  eigenvalue  k  defines  the  radial  variation  of  the  local  stress  field,  px . 
Solutions  of  k  covering  a  wide  range  of  fiber-matrix  combinations  (  0.0001<  (matrix 
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elastic  modulus  /fiber  modulus  )  <  10000)  are  shown  in  Figure  2.  These  results  are  for 
Poisson’s  ratios  of  vm=  vf  =0.3.  The  solution  demonstrates  that  X  is  real-valued  when 
the  fiber  is  stiffer  than  the  matrix.  For  some  fiber-matrix  combinations,  there  arejbd'Wo 
solutions  of  X  in  the  admissible  region  of  - 1  <  X  <  0.  We  denote  these  two  solutions  as 
Xl  and  X2,  (M  <  X2).  In  these  cases  there  are  two  singular  stress  terms  in  the  local 
analysis,  giving  rise  to  a  two-term  solution  : 

&°c*(p,d)  =  K{px'd\a-,6)+Kjx'd\k\6)+  ...  (10) 

Of  these  two  singular  terms,  the  most  negative-valued  X  (i.e.  X.j)  gives  rise  to  the  most 
singular  stress  distribution  at  O. 


Figure  3  depicts  that,  in  uncommon  ceramic  matrix  composite  materials  in  which 
the  matrix  is  very  much  stiffer  than  the  fiber,  X  can  be  complex-valued.  Since  complex 
roots  come  in  conjugate  pairs,  there  are  at  least  two  solutions,  Ai  and  X  \  ,  in  the 
admissible  range  -1  <  Re(A)  <  0.  For  complex-valued  X,  the  rank  of  the  system  of 
governing  equations  (5)  is  such  that  for  each  X  the  coefficients  in  X  can  be  determined  to 
within  one  arbitrarily  complex-valued  constant.  The  local  singular  stress  field 
corresponding  to  Xi  and  X  \  have  the  form: 


tf*(p  =  <C'(p  .e)-io%°‘{p  ,0) 

=  Kxpk'a^  (A, ;  0)  +  Kxp1'  of  (A, ;  0) 

(11) 

T^(p,0)  =  O,0)  +  <;>,0) 

=  K,pk'  ocg'  (A, ;  0)  +  Kxpk'  of'  (A, ;  0) 


The  singular  stress  field  of  (9)  or  (10)  is  a  locally  dominant  portion  of  the 
complete  stress  distribution.  The  complete  stress  field  (keeping  all  higher  order  terms  in 
(1))  has  the  form: 
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<^j(P^)  =  X  KJPlm<%&m'0>  (when  A  is  real)  (12) 

ms  1 

and 


r e(p  ,6)  =  <JM(p  ,6)-iope(p  ,0) 

=  2  am  ■,o)+KnpK  <?*-  dm ;  e) 

m  =  l 

zp(p  ,0)  =  opp(p  ,d)  +  iope(p  , 6) 

=  S  'f-P*'  oj"  u.;fl)  ■ +  kJ-  a;-  (!.;«) 


(when  X  is  complex)  (13) 


The  solution  of  the  local  analysis  is  valid  for  all  externally  applied  loading  and 
remote  geometric  constraints.  Different  externally  applied  loading  and  remote  geometric 
constraints  will  only  cause  difference  in  the  magnitude  of  the  scaling  factor  K. 


The  complex  potential  asymptotic  expansion  method  has  been  used  by  several 
researchers  to  study  stress  singularities  in  wedge  geometry  (Dempsey  and  Sinclair  (1981), 
Carpenter  (1985),  and  Iancu,  Fett  and  Munz  (1990)).  These  studies  focused  on 
determining  the  exponent  of  the  stress  singularity,  X  ,  at  the  wedge  comer.  The 
amplitude  factor  K,  and  hence  the  actual  magnitude  of  the  local  field  near  the  wedge 
comer,  were  not  solved.  There  are  other  good  methods  to  analyze  bonded  wedges,  see 
for  example  Bogy  (1971)  and  Hein  and  Erdogan  (1971).  Many  of  these  methods  are 
particularly  suitable  for  certain  class  of  global  loading  and  boundary  conditions,  and  are 
not  easily  extendible  to  determine  the  amplitude  factor  K  for  arbitrary  mechanical  and 
thermal  loadings.  In  the  local-global  matching  method,  the  amplitude  factor  K  for 
arbitrary  mechanical  and  thermal  loadings  for  arbitrary  global  boundary  constraints  can 
be  easily  determined. 

3  Global  Analysis  and  Local-global  Matching  Method 

In  the  global  analysis,  finite  element  method  (FEM)  is  used  to  determine  the 
complete  full-field  stress  distribution.  Then,  the  local-global  matching  method  is  used  to 
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link  the  full-filed  stress  distribution  to  the  local  field  to  determine  the  amplitude  factor  K. 

The  general  methodology  is  illustrated  here  with  one-term  matching  (matching  of  the 
one-term  local  field  (9),  at]  =  Klpx,dl  (X;8),  with  finite-element  results)  to  determine 

K|.  In  Section  4,  the  matching  method  is  extended  to  two-term  matching  to  determine 
both  Ki  and  K2  for  real-valued  Xj  and  X2.  In  Section  5,  the  method  is  generalized  to 
determine  the  complex-valued  scaling  factors  for  problems  with  complex-valued  Xs. 

An  axisymmetric  finite  element  model  with  quadratic  serendipity  (8-noded) 
quadrilateral  elements  were  used  to  compute  the  global  full-field  solution  (Figure  3). 
Near  the  fiber-protrusion  point  O,  elements  were  arranged  in  concentric  rings  with  the 

=  10“8(l  +  100i|n_11)  where  n  denotes  the  n-th  ring  of 

elements.  The  first  ring  of  elements  encircling  the  fiber-protrusion  point  had  an  outer 
radius  p/RQ.=  1.0  xlO  .  It  consisted  of  quadratic  triangular  elements  degenerated  from 

the  quadrilateral  elements  which  formed  the  rest  of  the  mesh.  No  singularity  shape 
function  or  special  implementation  for  stress  singularity  was  incorporated  into  the  finite 
elements.  The  0  -size  of  the  concentric  elements  was  22.5°.  The  mesh  consisted  of  389 
elements.  The  finite  element  computations  were  performed  using  the  ABAQUS  program 
(Hibbitt,  Karlsson  and  Sorensen,  Inc.  (1991)). 

Within  the  spatial  domain  of  the  finite  element  mesh,  we  searched  for  a  sub- 
region  where  the  angular  variation  of  the  axisymmetric  FEM-computed  complete  stress 
field  ,  a  rcM(0),  has  identical  0-dependence  as  that  of  the  plane-strain  singular  field, 
6^(0),  from  local  analysis.  If  the  angular  variation  of  FEM-computed  stresses  matches 

perfectly  with  angular  variation  of  the  singular  local  field  (9),  then  the  FEM-computed 
full-field  solution  has  captured  the  eigenvector  X  and  the  associated  eigenvalue  Xj  of 

the  local  singular  field.  Consequently,  inside  the  sub-region,  the  FEM-computed 
stresses  cr,^M  behave  with  the  spatial  structure  of  the  dominant  local  singular  field  (9): 

ofM  =  (14) 


ring  boundaries  given  by:  — 
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Therefore,  in  this  subregion  one  may  choose  any  stress  component,  say,<7w,  at  one 
point  at  any  angle,  say  0*.  to  determine  Kj , : 


*,= 


(p')V»(e’) 


(15) 


The  procedure  is  illustrated  here  through  an  example  composite  material.  The 
properties  used  correspond  to  those  of  Silicon  Carbide  (AVCO-SCS6  SiC)  fiber  and 
Borosilicate  matrix.  The  material  properties  and  the  eigenvalue  X  are  given  in  Table  1. 
The  local-global  matching  was  conducted  at  a  radius  p/R0.=  5.23  xlO'^.  For  an 

externally  applied  load  of  1  MPa  pulling  at  the  protruded  fiber  end,  the  value  of  Ki 
determined  was  0.1 145  MPa-mm"\ 

In  our  work,  we  determined  K[  from  cr^M  values  at  the  integration  point 
closest  to  the  fiber-matrix  interface  (at  0*  =  1.26  °).  If  the  finite  element  solution  has 
indeed  captured  the  characteristics  of  the  local  field  given  by  (9),  then  the  Ki  thus 
determined  from  one  stress  component  at  one  0  value,  will  make  all  stress  components 
of  the  finite-element  results  at  all  angles  0  to  automatically  match  with  the  singular  field 

of  (9).  Figure  4  compares  the  angular  variation  of  the  finite-element  solution  with 
0,^(0)  of  the  local  singularity  solution.  The  discrete  symbols  represent  finite  element 

results  and  solid  lines  denote  the  one-term  local  singularity  (9).  That  the  finite  element 
results  of  all  stress  components  at  all  0  values  indeed  match  the  singularity  solution, 
verifies  the  theory  and  assures  the  accuracy  of  Ki  so  determined.  This  local- 

global  matching  method  by  matching  the  entire  angular  variation  of  all  stress  components 

is  tantamount  to  matching  both  the  eigenvector  and  eigenvalue  of  the  local  field. 

For  each  of  the  three  stress  components  (<*pp>  099  and  Cpg),  we  define  a 

deviation  parameter  Djj(0)  to  measure  the  deviation  of  the  angular  variation  of  FEM- 
stress  from  that  of  the  one-term  local  field,  0^(0): 


-«!<•) 

kap) 


(16) 
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At  p/J\,=  5.23  xl(T^,  the  data  in  Figure  4  suggest  Djj(0)  is  close  to  zero  for  all 
components  of  stresses  at  all  0.  Similarly  results  were  obtained  for  the  several  other 
composite  systems  studied. 

Since  the  stress  field  has  a  singularity  at  O,  but  the  finite  elements  surrounding  O 
do  not  have  special  shape  function  to  simulate  unbounded  stress,  the  matching  is 
expected  to  be  less  than  satisfactory  within  the  very  first  ring  of  elements  surrounding  O. 
However,  starting  from  the  second  ring  of  elements,  results  show  good  matching 
characteristics  typical  of  Figure  4.  It  is  expected  that  as  ^  increases  beyond  a  certain 


extent,  the  one-term  singular  solution  (9)  deviate  more  from  the  complete  full-field 
stresses,  resulting  in  less  spectacular  angular  matching  at  sampling  points  at  large 

We  can  quantify  this  by  defining  a  mismatch  parameter,  M(0),  as  the  root-mean-square 
of  the  deviation  of  the  angular  variation  of  the  FEM-computed  stress  field  from  the 
angular  variation  of  the  one-term  local  field  (9): 


M(0)  = 


(17) 


In  (17),  the  following  convention  is  used:  <^=Opp,  and  °3=op0  ■  M(0)  is 

the  root  mean  square  of  the  three  D,j(0)  components. 

The  value  of  M(0)  at  different  radii  p  is  shown  in  Figure  5.  Clearly,  the  angular 

mismatch  depends  not  only  upon  the  radial  coordinate  p  but  also  depends  upon  the 

angular  coordinate  0  .  The  growth  of  the  mismatch  M(0)  with  5—  is  due  to  the 

Ko 

increasing  contribution  of  the  higher  order  terms,  (>-2.  X3,...),  to  the  complete  stress  field 

(12)  at  larger  . 

Ko 

To  measure  only  the  dependence  on  ,  we  construct  an  error  parameter  E] 

Ko 


from  M(0): 
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1  v 
e  =±y 
Nfi 


iy 

)|3* 


f3<£U<9,> 

Kpx  p 


(18) 


v '  J 

El  is  a  very  non-forgiving  measure  of  mismatch  error.  Only  perfect  matching  in  all 

stress  components  at  all  angular  positions  will  lead  to  a  zero  value  for  Ej.  Figure  6 

JL  _ m-8  _  A 

> 

example  composite.  This  plot  reveals  the  region  10"8  <  <  10'^  as  a  reasonable 


depicts  the  variation  of  E|  as  a  function  of  in  the  range  10'^  <  <  1  for  the 


domain  to  perform  one-term  matching  for  this  particular  composite  system  under  a  fiber 
pulling  load  of  1  MPa.  The  deviation  parameter,  mismatch  parameter  and  error 
parameter  (Djj(9),  M(0)  and  Ej  respectively)  all  measure  mismatch  but  give  different 

amount  of  detail.  They  are  useful  monitoring  parameters  to  quantify  the  quality  of  local- 
global  matching,  and  give  assurance  to  the  accuracy  of  Kj.  They  also  locate  the 

appropriate  Jr-  domain  within  which  the  one-term  local  field  is  dominant.  With  Kj 
Ko 


accurately  determined,  the  dominant  stress  distribution  (9)  at  the  fiber  protrusion  point  is 
completely  solved. 

The  local-global  matching  procedure,  and  the  associated  monitoring  parameters, 
are  also  applicable  to  problems  with  material  nonlinearity  (Lau,  Argon  and  McClintock 
(1983),  and  Tan,  Rahman  and  Lau  (1988)).  In  material  nonlinear  problems,  the  linear 
local  analysis  described  in  Section  2  has  to  be  replaced  by  asymptotic  solutions  for 
nonlinear  materials. 


4  Two-term  Matching 

The  example  composite  system  of  SiC(AVCO  SCS6)-Borosilicate  has  two 
singular  terms,  and  X.2-  In  regions  away  from  the  fiber  protrusion  point,  the 

dominance  of  the  one-term  field  (9)  reduces  and  higher  order  terms  start  to  gain 
influence.  The  two-term  local  field  (10)  has  a  larger  dominant  region  over  the  one-term 


14 


Pochiraju.  Lau  and  Wang 


field  (9).  Although  at  p/R0  =  5.23  xlO'^,  all  components  of  Djj(9)  are  close  to  zero 
(Figure  4),  at  larger  p/R0  the  Dy(0)  plot  reveals  an  interesting  feature.  In  Figure  7,  the 
discrete  symbols  denote  the  three  components  of  Djj(0)  at  p/Rg^  1.57  xlO'^.  The  lines 
denote  d ,/(0),  the  angular  variation  of  the  second  singular  term  (X2=-0.0092).  This 
figure  indicates  that  at  p/RQ.=  1.57  xlO’^the  complete  stresses  field  (as  given  by  the 

FEM-computed  stresses)  has  captured  more  features  of  the  two-term  singular  field,  (10), 
than  just  the  one-term  singular  field,  (9).  In  fact,  if  the  FEM-computed  stresses  exactly 
match  the  two-term  singular  field,  according  to  (10)  and  (16)  D[j(0)  should  match  with 
d^{9)  within  an  arbitrary  scaling  constant.  This  observation  motivates  a  general  two- 

term  local-global  matching. 

Two  different  methods  have  been  used  to  obtain  both  Ki  and  K2  with  equal 
success.  The  first  method  uses  directly  the  two-term  solution  in  the  local  analysis  (10). 
Two  stress  components  at  one  spatial  point  or  one  stress  component  at  two  different 
spatial  points  are  matched  and  a  system  of  two  linear  equations  is  solved  to  obtain  the 
amplitude  factors  Ki  and  K2-  In  this  study,  two  different  spatial  points  along  the 
interface  (same  0  but  different  p)  are  chosen  and  the  normal  stresses  are  extracted  from 
the  global  finite  element  analysis.  Let  the  two  points  be  at  radii  piand  p2,  then 

<^m(a.o)=*.aa,+*2AAi 

^(Pi,o)=k,p^+k,p21’'  (19) 

These  two  equations  are  solved  simultaneously  to  give  K  i  and  K2- 

The  discrete  symbols  in  Figure  8  denote  the  normalized  FEM-computed  Oqq  at 

0  =2.6°  (close  to  the  fiber-matrix  interface)  at  varies  p/Ro.  The  finite  element  results 
theoretically  represent  the  complete  full-field  stress  distribution.  The  radial  variation  of 
the  FEM-stresses  are  compared  with  the  one-term  singular  field  shown  in  dashed  line, 
and  the  two-term  singular  field  shown  in  solid  line.  The  one-term  field  is  calculated 
using  (9)  with  Ki  obtained  from  one-term  matching.  The  two-term  field  is  calculated 
using  (10)  with  K(  and  K2  obtained  from  two-term  matching.  In  regions  close  to  the 
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fiber  protrusion  point,  the  one-term  field  provides  a  good  matching  with  the  FEM- 
computed  stresses  (see  also  Figure  6  at  p/Ro=2.53  x  10'^).  At  distances  away  from  the 
fiber  protrusion  point  (say,  at  p/Ro  =  10  ),  using  (9)  with  the  one-term  matching  does 
not  acceptable  approximation  to  the  complete  stress  solution  (represented  by  FEM 
results).  Figure  8  shows  that  (10)  with  the  two-term  matching  results  give  reasonable 
approxmation  to  the  complete  stress  distribution  up  to  p/Ro  =  10' ^ 

In  Figure  9,  finite  element  results  are  compared  with  the  one-term  singular  field 
(9)  and  two-term  singular  field  (10).  The  discrete  symbols  are  the  FEM-computed 
angular  variation  of  stresses  at  p*/Ro=3.33  x  10^.  The  dashed  line  denote  the  angular 
variation  a1  (A,;0)  of  the  one-term  field  (9).  The  solid  lines  denote  Sjj(0)  given  in 

(20).  Components  of  Sy(9)  represent  the  ngular  variation  of  stresses  of  the  two-term 
field  (10)  at  the  radius  p  =  p*. 


S„(0)  = 


' ATjfp*)*1  dj,(0)  +  K2(p)X:  crv(0)' 

«rv.o) 


(20) 


The  two-term  matching  procedure  was  applied  to  the  Silicon  Carbide-Borosilicate 
composite,  for  two  different  types  applied  loading:  that  of  fiber  pulling  with  1  MPa,  and 
that  of  thermal  loading  by  a  change  of  temperature  of  DT=  -1°C.  Table  2  lists  solutions 
of  Kj  and  K2  determined  by  the  two-term  matching  procedure.  Also  shown  are  results 
for  another  composite  system  (Nicalon  fiber  -Silicon  Carbide  matrix  with  Ef  =  179  GPa, 
Em  =379  GPa,  vf  =  0.2  and  vm  =  0.2).  The  external  loading  for  this  latter  case  is  IMPa 
of  applied  fiber-pulling  traction. 

With  K]  and  K2  so  determined,  the  locally  dominant  singular  field  is  completely 

known  (10).  The  relative  dominance  of  the  first  and  second  term  in  the  two-term  field 
depends  upon  the  values  of  A|  and  and  the  applied  loading.  The  two  parameters  that 

govern  this  relative  dominance  are:  the  difference  between  the  two  eigenvalues,  (Aj-  Aj), 

j^2 

and  the  ratio  of  the  two  amplitude  factors  .  The  first  parameter  can  be  solved  by 
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local  analysis  alone,  but  the  latter  parameter  is  governed  by  the  applied  loading  and  must 
be  determined  by  two-term  local-global  matching. 

The  second  method  for  two-term  matching  is  based  on  (10)  and  (16).  In  the 
region  in  the  finite  element  mesh  where  the  two-term  singular  field  (10)  is  dominant. 


<r,ravVe)  -  k.pV.W  ■  kj>1- aj(e) 


(21) 


and 


0„(«)  *  § HpT'1'  «5(«) 

K\p  ) 


(22) 


In  this  procedure,  Ki  is  first  determined  as  described  in  the  Section  3  from  a  sub-region 
where  the  one-terms  singular  field  dominates.  Then,  in  a  further  region  in  which  the 
two-term  singular  field  dominates,  K2  can  be  determined  from  angular  variation 
matching  d~(6)  using  (21)  directly.  Or  the  ratio  K2/K1  can  be  determined  through 

angular  variation  matching  using  (22).  Using  (21)  or  (22),  K2  can  be  obtained  with  the 
same  procedure  for  one-term  matching  of  Section  3. 


5  Matching  for  Complex-valued  Stress  Singularity 

The  methodology  developed  in  this  paper  is  applicable  to  the  cases  of  complex 
stress  singularities  with  equal  merit.  To  illustrate  this  fact,  a  hypothetical  composite  with 
the  matrix  modulus  fifteen  times  greater  the  fiber  modulus  is  considered.  The  material 
properties  and  the  two  complex-valued  eignevalues ,  and  X  \  ,  are  listed  in  Table  3. 

The  complex-valued  scaling  factor,  K,  for  problems  with  complex-valued 
singularity  exponent  X{  and  X  {,  is  defined  by  the  following  equation. 


Ta(p\0)  =  c™(p\0)-icr»(p,0) 

=  K(p)X 


(23) 


or 


K  = 


r»(p*.Q) 

(P? 


(24) 
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The  error  parameter  Ej  (18)  can  be  generalized  to  the  case  of  complex-valued  X 
as  follow. 


in  which 

|/fpA|  =  1/(Arpi  (26) 

The  finite  element  mesh  used  for  this  case  of  complex-valued  stress  singularity  was  the 

same  as  described  in  Section  3.  The  finite-element  computed  interfacial  normal  and 

shear  stresses  are  used  in  the  local-global  matching.  In  this  example,  the  interfacial 

normal  and  shear  stresses  were  determined  at  =3.884  x  10'^.  (Note  that  once  this 

Ko 

radius  is  chosen,  the  complex  coefficients  in  the  local  analysis  have  to  be  normalized 

such  that  the  angular  variation  function  at  the  interface  at  this  radius  is  unity,  (23).  To 
achieve  this,  the  complex-coefficient  Af  was  set  to  a  value  of  (-0.013098  +  i  0.018136) 

in  the  local  analysis. 

In  the  local-global  matching,  the  complex-valued  stress  amplitude  factor  was 

computed  from  (24),  giving  K  =  (- 0.0887  -  i  0.0183563)  corresponding  to  the  stress 

singularity  exponent  Xi  =(-0.456  +  i  0.07291).  The  angular  variation  matching  was 

found  to  be  very  satisfactory.  Figure  10  shows  that,  with  the  complex-valued  K  so 

determined  from  one  0*  value,  the  singular  field  (denoted  by  lines)  from  local  analysis 

agreed  well  with  the  FEM-computed  stresses  (discrete  symbols)  at  all  values  of  6  for 

all  components  of  stresses.  Figure  1 1  shows  the  variation  of  the  error  Ej  (24)  with  Jr- . 

Ko 

Comparing  Figure  1 1  with  Figure  6  (the  one-term  matching  in  a  local  field  where  there 
two  negative  Xs),  it  seems  that  Ei  for  the  complex-valued  X  is  comparatively  less  at  large 

p 

s—  .  The  main  reason  for  this  is  that  the  only  singular  terms  in  the  complex-valued  case 
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are  due  to  Xi  and  X).  Since  Xi  and  X  i  are  complex  conjugates  with  identical  real 

parts,  this  complex-valued  matching  has  incorporated  all  the  singular  terms  in  the  local- 

global  matching.  On  the  contrary,  the  one-term  matching  results  of  Figure  6  has  not 

considered  the  other  admissible  singular  term  ,  X2.  Therefore,  at  large  5—  the  neglected 

Ko 

X2  term  contributes  to  Ej  in  Figure  6. 

6  Comparison  with  Purely  Finite  Element  Estimation 

Without  conducting  local  analysis  to  determine  Xand  a y(0),  but  assuming  the 
local  stresses  is  a  one-term  singular  field  of  the  form:  <7i;(p,0)  =  K  pxdtI{d)  where  X  is 
real-valued,  one  may  use  finite  element  results  alone  to  estimate  X  and  K.  The  Finite 
element  estimation  techniques  generally  involve  finding  the  slope  of  the  radial  variation 
of  the  stress  components  on  a  log-log  scale.  Figures  12(a)  and  12  (b)  show  such  log-log 
plots  to  estimate  X  using  the  finite  element  results  from  our  global  analysis.  From 
Figures  12(a)  and  (b),  it  is  clear  that  the  value  of  X  so  estimated  from  purely  finite 
element  results  varies  in  magnitude  depending  on  which  stress  component  is  used  and 
which  angle  0  is  used. 

The  stress  amplitude  factor  K  can  also  be  estimated  using  purely  finite  element 
results.  In  such  estimation  of  K  one  needs  First  the  value  of  X.  Since  the  X  estimated 
from  finite  element  results  has  a  range  of  values  depending  on  the  stress  component  used 
or  the  0  used,  this  uncertainty  is  passed  on  to  the  value  of  K  when  K  is  estimated.  There 
is  no  good  means  to  quantify  the  error  in  estimating  K  and  X  using  purely  finite  element 
results.  Without  solving  the  local  singularity  a  priori,  there  is  little  assurance  on  the 
accuracy  of  K  estimated  from  the  Finite  element  results. 

On  the  other  hand,  a  prior  knowledge  of  the  singularity  exponent  X  helps  guide  to 
determine  K  accurately.  For  crack  problems,  knowing  that  X  =  -0.5  helps  guide  the 
use  of  Finite  element  results  to  solve  for  stress  intensity  factor.  For  the  free  edge 
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interface  in  bimaterials,  Munz  and  Yang  (1992)  first  solved  X  by  Bogy’s  (1971)  approach 
to  guide  subsequent  use  of  finite  element  results  to  determine  K  accurately.  For 
problems  where  X  is  not  known,  it  is  recommended  to  conduct  first  a  local  analysis  before 
using  finite  element  results.  In  the  local  analysis,  other  than  calculating  just  X  it  is  worth 
a  researcher's  time  to  determine  also  the  angular  variation  of  stresses.  With  both  the  X 
and  angular  variation  known,  in  subsequent  local-global  matching  the  researcher  can  seek 
to  match  both  the  eigenvalue  and  eigenvector  to  ensure  more  reliable  results. 

Application  of  purely  finite  element  results  in  log-log  plots  to  estimate  X  and  K 
does  not  work  well  to  solve  for  the  two-term  singular  field  (10)  should  it  exists. 
Although  the  example  composite  system  Silcon  Carbide-Borosilicate  has  two  singular 
terms  X\  and  X2  (Table  2),  the  finite  element  results  does  not  show  any  change  in  slope 

(Figures  12(a)  and  (b)).  Consequently,  without  local  analysis,  just  the  finite  element 
results  would  not  have  enough  resolution  to  solve  for  both  Xj  and  X2  ,  and  both  Kj  and 

K2. 

When  the  stress  field  has  complex-valued  X,  then  finite  element  estimation 

techniques  based  on  radial  variation  of  stresses  (Figures  12(a)  and  12(b))  would  not 
work.  With  complex-valued  X,  a  typical  stress  component ,  say  090,  is  of  the  form 

age  =  Re(re)  =  Re[tfpAFtf(cos(£logp),sin(£logp),A,0)j,  (where,  e  =  Im{A] ) 

which  shows  stress  oscillation  with  respect  to  p.  On  the  other  hand,  the  local-global 
matching  method  presented  in  previous  sections,  is  applicable  to  stress  fields  with 
complex-valued  X. 

A  relative  advantage  of  the  local-global  matching  method  is  that  the  monitoring 
parameters  Dy(0),  M(0)  and  Ej  can  be  used  to  measure  the  quality  of  matching.  Use  of 

these  quality  control  parameters  can  ensures  an  accurate  determination  of  K  by  verifying 
that  the  finite-element  computed  stresses  used  in  matching  indeed  has  captured  the  full 
spatial  characteristics  of  the  asymptotic  singular  solution  in  the  space  (domain)  of 
matching. 
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7  Concluding  Remarks 

The  local  analysis  presented  can  determine  the  detail  spatial  structure  of  singular 
stress  field  at  the  fiber  protrusion  point.  For  the  problems  considered  here,  both  one-term 
singular  field  and  two-term  singular  field  are  possible.  In  some  composite  materials,  the 
singularity  exponent  can  be  complex-valued. 

The  local-global  matching  method  can  determine  the  amplitude  factors  accurately 
for  the  one-term  field  and  two-term  field  with  real-valued  singularity  exponent.  It  is 
based  on  angular  variation  matching  which  is  tantamount  to  matching  both  the  eigenvalue 
and  eigenvector  of  the  asymptotic  field.  It  is  recommended  that  in  local  analysis,  not 
just  the  singularity  exponent  X  but  also  the  angular  variation  of  stresses  should  be 
determined  so  one  can  match  the  eigen-pair  instead  of  just  matching  the  eigenvalue  in 
subsequent  local-global  matching. 

The  monitoring  parameters  Djj(0),  M(0)  and  Ej  can  provide  quantifiable 
measures  of  the  quality  of  local-global  matching. 

In  its  generalized  form,  the  local-global  matching  is  equally  applicable  to  cases 
with  complex-valued  singularity  exponent  to  determine  complex-valued  amplitude  factor. 

The  local-global  matching  method  can  be  applied  for  other  interfacial  conditions 
such  as  frictional  sliding  and  interfacial  crack. 

With  minor  modification,  the  local-global  matching  method  is  also  applicable  to 
material  nonlinearity  problems. 
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Table  1:  Material  properties  and  singularity  exponents,  Xi  and  X2,  for  the  Silicon 
Carbide-Bonosilicate  system. 


IB 

IB 

vm 

Vf 

am-af 

xlO'5 

Fiber  Radius 
(Ro) 

A-l 

>■2 

63  GPa 

410  GPa 

0.2 

0.2 

7.2  /°C 

71.5  pm 

-0.313344 

-0.0092246 

I 

I 


Table  2:  Amplitude  factors  for  the  first  and  second  singular  terms  for  three  example 
problems. 


Material/ 
(Loading  Type) 

Xl 

>-2 

m 

m 

Sic-Borosilicate 
(Pulling- 1  MPa) 

-0.313344 

-0.00922466 

0.1153 

-0.2998916 

Sic-Borosilicate 
(DT=  -  IOC) 

-0.313344 

-0.00922466 

0.0333911 

-0.12094 

Ni-  Sic 

(Pulling- 1  MPa) 

-0.4914585 

-0.1747356 

0.01458 

-0.05103 

Table  3:  Material  properties  of  the  hypothetical  composite  system  having  complex-valued 
singularity  exponents,  A.1  and  X 1 


■1 

■1 

vm 

Vf 

WBSSm 

Fiber  kadius 
Ro 

15  GPa 

1  GPa 

0.2 

0.2 

1  /°c 

71.5  p.m 

(-0.456  +  i 

0.07291) 

(-  0.456  -  i 

0.07291) 
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Legends  to  Figures 

Figure:  1.  Schematics  of  the  single  fiber  pullout  model  showing  local  and  global  co¬ 
ordinate  systems. 

Figure  2:  Variation  of  the  exponent  of  stress  singularity  X  with  material  properties 
Figure:  3.  Schematics  of  the  finite  element  mesh 

Figure  4:  Comparison  of  angular  variation  of  FEM-computed  stresses  (discrete  symbols) 

O 

with  one-term  singular  field  (lines)  at  p/Ro=  5.3x10'°  for  an  example 
composite  Silcon  Carbide-Borosilicate  (SiC-Borosilicate). 

Figure  5:  Variation  of  the  mismatch  parameter  M(0)  with  0  in  one-term  matching 

Figure  6:  Variation  of  the  error  parameter  E]with  p/R0  in  one-term  matching. 

-  i  -4 

Figure  7:  Comparison  between  Dy(0)  and at  p/Ro=l. 57x10  . 

Figure  8:  Comparison  of  FEM-computed  full-field  stress  (discrete  symbols)  with  one- 

term  singular  field  (dashed  line)  and  two-term  singular  field  (solid  line)  for  the 
SiC-Borosilicate  example  composite. 

Figure  9:  Comparison  of  FEM-computed  angular  variation  with  dr1  (A, ;  6)  and  Sjj(0)  at 

'l 

r/R0=3. 33x10' J  for  the  SiC-Borosilicate  example  composite. 

Figure  10:  Angular  variation  matching  for  composite  with  complex-valued  stress 
singularity  exponent,  t/Rq=  5.28x10'^. 

Figure  11:  Variation  of  the  error  parameter  E^with  p/Ro  for  composite  with  complex¬ 
valued  stress  singularity  exponent. 

Figure  12  (a):  Estimation  of  the  exponent  of  singularity  X  based  radial  variation  of  finite 
element  results  of  two  stress  components. 

Figure  12  (b):  Estimation  of  the  exponent  of  singularity  X  based  radial  variation  of  finite 
element  results  of  am  at  three  different  angles. 
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Figure  2:  Variation  of  the  exponent  of  stress  singularity  X  with  material  properties 
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Figure:  3.  Schematics  of  the  finite  element  mesh 


Pochiraju,  Lau  and  Wang 


0  (degrees) 


Figure  4:  Comparison  of  angular  variation  of  FEM-computed  stresses  (discrete  symbols) 
with  one-term  singular  field  (lines)  at  p/Ro=  5.3x10'°  for  an  example 
composite  Silcon  Carbide-Borosilicate  (SiC-Borosilicate). 
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Figure  5:  Variation  of  the  mismatch  parameter  M(0)  with  0  in  one-term  matching 
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Figure  7:  Comparison  between  Djj(0)  and  dl}2(6)  at  p/Ro=  1.57x10 
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Figure  8:  Comparison  of  FEM-computed  full-field  stress  (discrete  symbols)  with  one- 
term  singular  field  (dashed  line)  and  two-term  singular  field  (solid  line)  for  the 
SiC-Borosilicate  example  composite. 
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Figure  9:  Comparison  of  FEM-computed  angular  variation  with  (A,;0)  and  Sjj(0)  at 
r/R0=3. 33x10'^  for  the  SiC-Borosilicate  example  composite. 
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Figure  11:  Variation  of  the  error  parameter  E  \  with  p/Ro  for  composite  with  complex¬ 
valued  stress  singularity  exponent. 
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Figure  12  (a):  Estimation  of  the  exponent  of  singularity  X  based  on  radial  variation  of 

finite  element  results  of  two  stress  components. 
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Figure  12  (b):  Estimation  of  the  exponent  of  singularity  X.  based  on  radial  variation  of 
finite  element  results  of  <J69  at  three  different  angles. 
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fhe  coefficient  of  thermal  expansion  (CTE)  of  a 
1 40-pin  -diameter  SiC  monofilament  was  determined  to  be 
6.5  ±  0.5  x  1 in  the  temperature  range  25°-900°C. 
Heat  treatment  of  the  fibers  at  1400°C  for  90  min  resulted  in 
both  an  increase  in  the  intensities  of  the  SiC  peaks  and  a 
reduction  of  their  width,  indicating  some  grain  growth. 
Furthermore,  heat  treatment  of  the  fibers  for  2  h  at  1600°C 
in  vacuum  reduced  their  CTE  to  5.7  x  lO'^C"1.  The 
anomalously  high  CTE  and  its  decrease  upon  heat  treat¬ 
ment  were  attributed  to  the  presence  of  a  large  fraction  of 
the  atoms  in  the  boundaries  between  the  nanometer  3-SiC 
grains.  The  thermal  expansion  coefficients  of  a  series  of 
SiC'borosilicate  composites  were  measured  as  a  function  of 
fiber  volume  fraction  in  the  temperature  range  25°-500°C 
and  found  to  follow  the  rule  of  mixtures. 

I.  Introduction 

Is  ni.siiiMv;  ceramic  ceramic  composites  one  oi  (he 
more  important  considerations  is  the  thermal  compatibility 
between  the  reinforcing  phase  and  that  of  the  matrix.  For  tiber- 
reinforced  ceramic  matrix  composites  both  the  radial  and  the 
axial  thermal  expansion  mismatches  are  important.  If  the  fibers 
have  a  larger  coefficient  of  thermal  expansion  (CTE)  than  the 
matrix,  then,  upon  cooling  from  the  processing  temperature, 
lhe  matrix  w  ill  be  in  axial  compression.  Radially  the  fiber 
matrix  interface  will  be  in  tension  and  will  tend  to  debond 
Conversely,  if  the  CTE  of  the  fibers  is  lower  than  the  matrix,  the 
latter  w  ill  radially  clamp  dow  n  on  the  fibers  upon  cooling,  but 
the  matrix  will  experience  an  axial  residual  tensile  stress.  In 
cither  case,  degradation  in  the  mechanical  properties  of  the 
composites  can  result. 

One  of  the  more  promising  fibers  for  reinforcing  ceramic 
matrices  is  the  SCS-6  SiC  monofilament  (Textron  Specialty 
Materials.  Lowell.  MA).  In  addition  to  its  high  intrinsic 
strength  and  stiffness,  it  is  fairly  oxidation  resistant.  This  fiber 
is  manufactured  by  chemical  vapor  deposition  of  SiC  onto  a 
35-p.m  carbon  core.  A  micrograph  of  the  cross-section  of  the 
fiber  is  shown  in  Fig.  I .  The  SiC  grains,  composed  mainly  of 
[i-SiC,  are  heavily  faulted,  radially  oriented,  and  columnar 
with  the  octahedral  planes  preferentially  onented  to  the  deposi¬ 
tion  surface.1  The  graphitic  carbon  coating  on  the  fiber  surface 
is  not  pure,  but  has  been  shown  to  have  SiC  particles  embedded 
in  it.1 

The  radial  and  axial  thermal  expansion  coefficients  of  SCS-6 
monofilaments  have  been  measured  by  various  investigators'  “ 
and  the  results  are  summarized  in  Table  I.  The  values  for  the 
axial  CTE's  range  from  4  x  IO',0C"'to5  5  x  I0'"5C  '1 
Although  part  of  the  scatter  can  be  attributed  to  the  different 
temperature  ranges  investigated,  the  range  is  quite  wide.  As 
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Fig-  l-  Micrograph  ot  SiC  monofilament. 


will  be  discussed  in  more  dctai I  below,  another  possible  source 
of  variability  may  be  differences  in  the  amorphous  phase  con¬ 
tent  between  different  manufactured  fiber  lots. 

The  reported  values  for  the  radial  CTE's  also  vary  (see 
Table  1)  As  far  as  can  be  traced,  the  only  direct  measurement 
was  conducted  by  Goettler  and  Faber."  Other  values  reported' T 11 
were  interred  from  stress  bifringenee  measurements  conducted 
at  room  temperature  as  (he  thermal  expansion  of  the  matrices 
were  v  aried  In  one  study.  ‘  the  decrease  in  intertacial  shear 
stress  required  to  pull  out  single  fibers  from  a  borosilicate 
matrix  as  the  temperature  increased  was  taken  as  evidence 
that  the  radial  CTE  of  the  fibers  was  less  than  the  matrix,  i.e.. 
<3.2  ■  10  "°C  1 


Table  i.  Summary  of  Reported  Values  of 
CTE  of  SCS-6  Fibers  and  Composites 


u  mp 

Rjnirc 
.  Ci 

Lit: 

•  urc 

Commem-s 

Rcl 

Axial 

400-1350 

5  2 

Taken  from  Fig  4  in  Ref.  3 

3 

6(H)- 1200 

5  0 

Taken  from  Fig.  5  in  Ret.  3 
Pronounced  shrinkage 
observed  at  I360'C 

8(H)- I3»H) 

5  5 

Taken  from  Fig.  7  in  Ret  4 

4 

25-5<H) 

3.99 

Dilatometry  (data  was 
not  shown). 

s 

25-vHH) 

6  5-05 

This  work 

R  ad i a  1 

2.63 

Dilatomeirv 

6 

4  0 

Based  on  stress  bitringence 
measurements 

s 

■  ;  2 

Based  on  stress  bitringence 
measurements 

•  3  2 

Based  on  stress  bifringenee 
measurements 

s 
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Measurement  of  relaxation  ot 
residual  stresses 
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Fig.  2.  Schematic  ol  sample  holder  used  tor  measuring  (he  axial  ther 
mal  expansion  ol  the  libers 


The  objective  of  this  work  was  to  measure  the  axial  thermal 
expansion  of  SCS-6  SiC  monofilament  and  composites  made 
using  these  fibers  in  an  attempt  to  better  quantify  and  under¬ 
stand  their  thermal  expansion  behavior. 

II.  Experimental  Procedure 

I  / )  Tiber  T her  mal  Expansion 

The  axial  thermal  expansion  ot  (he  SiC  monofilaments 
i  grade  SCS-b.  AYCO  Specialty  Materials.  Textron.  Inc  .  Low¬ 
ell.  MAI  was  measured  using  the  sample  holder  shown  sche¬ 
matically  m  Fig.  2.  To  align  (he  fibers  and  keep  them  from 
buckling.  H— 10  libers  ( 10-12  mm  long)  were  inserted  through 
holes  drilled  in  the  centers  of  ceramic  disks  machined  out  ot 
machinable  ceramics,  such  that  =2  mm  of  fibers  protruded  on 
either  side.  The  disks  in  turn  were  placed  in  a  dilatometer 
( Model  No  TDA-HI  with  an  alumina  holder,  Harrop  Indus¬ 
tries.  Columbus,  OH)  and  the  thermal  expansion  was  measured 
in  the  temperature  range  between  room  temperature  and  900°C 
The  fibers  were  tree  to  slide  through  the  ceramic  disks  at  all 
times  and  no  reaction  was  observed  between  the  fibers  and  the 
disks 

The  sample  temperature  was  measured  using  a  thermocouple 
placed  as  close  as  possible  to  the  fibers.  Heating  the  fibers  at 
temperatures  10(M)oC  resulted  in  significant  shrinkage.  This 
shrinkage  is  believed  to  be  due  to  buckling  and/or  creep  of  the 
libers  under  the  push-rod  loading. 

Preliminary  runs  indicated  that  good  thermal  contact 
between  the  fibers  and  the  dilatometer  fixture,  together  with 
slow  heating  rates,  gave  the  most  reproducible  results.  The 
heating  rale  was  5’C  min  To  ensure  good  thermal  contact,  as 
well  as  completely  shield  the  fibers  from  oxidation,  the  whole 
assembly  was  covered  with  a  graphite  powder  as  shown  in 
Pig.  2  Ar  gas  was  passed  over  the  samples  at  all  times.  All  the 
data  presented  here  were  obtained  with  the  C-cover  and  under  a 
flow  ing  atmosphere  of  Ar.  Nevertheless,  the  thermal  expansion 
at  temperatures  below  200CC  were  anomalously  low;  the  reason 
lor  this  is  not  entirely  clear  but  could  possibly  be  due  to  a  ther¬ 
mal  lag  between  the  thermocouple  and  the  fibers.  This  behavior 
is  not  uncommon  for  this  type  of  dilatometer. 

The  dilatometer  was  calibrated  by  measuring  the  thermal 
expansion  ol  fused  quartz,  a  borosilicatc  glass  (7740.  Coming 
Glass  Works.  Coming.  NY),  and  a  monolithic,  fully  dense  SiC 
cylinder.  In  all  cases  the  measured  values  agreed  with  the  pub¬ 
lished  results. 

(2)  Thermal  Expansion  of  Composites 

Unidirectional  reinforced  composites  were  fabricated  by  vac¬ 
uum  hot-pressing  borosilicate  glass  powders  and  the  SiC  mono¬ 
filaments  at  I200°C  for  10  min  The  details  of  the  fabrication 
procedure  can  be  found  elsewhere  10  The  fiber  volume  fraction 
of  each  specimen  was  individually  measured  by  determining 
the  total  fiber  area  in  the  specimen's  cross  section.  Bars.  10-20 
mm  long,  of  square  cross  section  (2  mm  x  2  mm)  were  used 
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Pig.  3.  Functional  dependence  ol  axial  thermal  expansion  ol  ihe  SiC 
monolilainenis  on  lemperature  in  the  25  -OIK)  range  The  lop  lour 
curves  arc  tor  tour  different  hber  batches.  The  middle  curve  (open 
squares)  is  lor  monolithic  dense  SiC  The  bottom  curve  is  for  7740 
borosilicate  glass 


tor  the  measurements  The  temperature  range  over  which  the 
thermal  expansion  could  be  measured  was  limited  on  the  high 
side  hv  the  dilatometric  softening  point.  (  =550°C)  of  the  bortv- 
vilicaie  glass  matrix  The  CTF.  ot  ihe  hot-pressed  glass  matrix 
was  also  measured  (see  Pig.  .'i  and  found  to  be  3. 1  »  10  "C. 
in  agreement  with  manufacturer's  data  sheet. 


III.  Results  and  Discussion 


( / )  Thermal  Expansion  of  Eibers 

The  temperature  dependence  of  the  fiber  displacements  in  the 
lemperature  range  25°-d(X)  'C  for  four  different  fiber  samples  is 
shown  in  Pig.  3.  The  corresponding  CTE  in  the  temperature 
range  200°-9<H)oC  is  6  5  *  0.5  x  10  'M,C ' 1  This  value  is 
higher  than  any  of  the  reported  values  to  date  (see  Table  1). 
Also  included  in  Fig  3  are  the  results  of  a  monolithic  block  of 
dense  SiC.  where  the  measured  CTE  is  calculated  to  be  4  5  X 
10  T  1  It  is  worth  noting  that  this  value  is  significantly  Iowa 
than  that  of  the  fibers. 

Treating  the  fiber  as  a  composite  w  ith  a  C-core  in  a  SiC 
matrix  ( V'.,,  =  K'T.  E  ,„  =  40GPa.  anda„„  =  10  x  10  *°C 
(Ref.  inand£s,<  =  440  GPa)  and  using  Eq.  1  (see  below),  the 
calculated  fiber  CTE  is  4.54  x  I0""°C  1  which  clearly  cannot 
account  for  the  CTE  value  measured  here.  The  thermal  expan¬ 


sion  amsotropvi:  of  SiC  is  also  not  sufficient  to  explain  the 
anomalously  high  a  of  the  fibers  _*~ 

X-ray  diffraction  ( XRD)  of  the  as-received  fibers  shows  then 
to  be  composed  of  fi-SiC.  However,  TEM  observations  earned 
out  by  Wawner  et  at."  14  have  shown  that  the  bulk  of  these  SC 
monofilaments  consisted  of  columnar  (3-SiC  subgrains  WJ?  <£ 
were  30-100  nm  in  size.  This  implies  that  a  large  fraction  of 
atoms  are  residing  in  the  grain  boundary  areas.  We  propose  “<■  < 
these  atoms,  which  can  be  considered  to  be  an  amorpho**®^ 
gram  boundary  ‘phase.”  are  responsible  tor  the  higher  CTfc^ 
the  fibers  as  compared  to  monolithic  SiC.  To  test  this  hy 
sis  a  number  of  fibers  were  covered  in  graphite  and  heat-) 
at  I400°C  under  flowing  Ar  for  90  min.  The  X-ray  po 
traction  of  the  as-received  and  heat  treated  samples.  * 

MgO  internal  standard,  is  shown  in  Fig  4.  « 

ments  were  repealed  on  three  different  batches  of 
similar  results.  Both  an  increase  in  the  intensity  o  -ad 
peaks  and  a  decrease  in  the  full  width  at  one-half  MM 
the  peaks  were  observed  as  a  result  of  the  heat  treatm®*^ 
results  are  clear  and  unambiguous  evidence  k5'.'**'-. . 
crystallinity  of  the  samples  and  the  occurrence  of  g  ^ 

From  the  X-ray  line  broadening,  using  ihe  ScherrW^j 
the  average  crystallite  size  was  calculated  to  be  2 
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Fig.  4.  X-ray  powder  diffraction  of  as-received  libers  (solid  dia¬ 
monds)  and  libers  (hat  were  heal  treated  at  I400X  lor  00  min  (open 
squares).  The  peak  on  the  right  is  that  of  the  M gO  standard 


as-received  fibers  and  50  nm  for  the  heat  treated  ones,  which  is 
in  good  agreement  with  the  values  reported  by  Wawner  et  al. 

An  implication  of  this  hypothesis  is.  that  the  amorphous  con¬ 
tent  of  the  fibers  influences  their  CTE's.  To  test  this  hypothesis, 
the  CTE  of  a  set  of  fibers  that  were  heat-treated  in  vacuum 
(=0. 1  Pa)  for  2  h  at  1500°C  was  measured.  The  CTE  dropped 
to  5.7  x  10  “C.  which  is  lower  than  the  value  for  the  as- 
received  fibers.  From  the  work  conducted  here,  what  remains 
unclear  is  the  distribution,  volume  fraction,  and  thermal  expan¬ 
sion  of  the  amorphous  phase.  Further  work  would  have  to  be 
conducted  to  explore  this  issue. 

There  is  indirect  evidence  in  the  literature  for  the  occurrence 
of  structural  changes  in  these  fibers  in  the  temperature  range 
900°-1400°C-  Based  on  creep  measurements  DiCarlo4  con¬ 
cluded  that  thermally  induced  microstructural  changes  that 
reduce  the  creep  rate  occurred  between  1200°  and  I400°C. 
Furthermore.  McHcnery  and  Tressler"  observed  a  strength 
degradation  in  the  900°- 1 ()00°C  range  and  speculated  that 
recrystallization  and/or  grain  growth  of  the  SiC  may  be  respon¬ 
sible.  However,  their  XRD  and  SEM  analysis  failed  to  identify 
a  specific  mechanism.  Crystallization  of  an  amorphous  grain¬ 
boundary  phase  and/or  grain  growth  could  explain  both 
observations. 

It  is  thus  possible,  that  during  manufacturing  variations  in 
processing  variables  could  have  resulted  in  variations  in  the 
grain  size  or  amorphous  content  of  the  fibers.  This  would  not 
only  explain  the  variability  in  the  CTE  values  measured 
between  various  investigators,  but  also  explain  the  anoma¬ 
lously  high  values  of  CTE  measured  here,  ft  would  also  explain 
the  variability  in  shrinkage  reported  from  lot  to  lot.' 


Fig.  5.  Effect  of  fiber  volume  fraction  on  the  axial  thermal  expansion 
coefficients  of  composites.  Dotted  line  is  a  plot  of  Eq  ( I )  (sec  text). 


(2)  Thermal  Expansion  of  Composites 

The  functional  dependence  of  the  axial  CTE  of  the  composite 
samples  on  fiber  volume  fraction  is  shown  in  Fig.  5,  where  each 
datum  point  is  the  result  of  at  least  two  or  more  separate  sam¬ 
ples  at  each  volume  fraction.  The  dotted  line  was  calculated 
according  to  the  following:1' 

__  a,E,C  +  «„,E„,V', 

u'""’p  ~  fc.v;  +  e,„v;„  1 1 1 

where  f  and  m  refer  to  the  fiber  and  matrix,  respectively,  a,. 
and  V,  are  the  CTF..  Young's  modulus,  and  volume  fraction, 
respectively.  Figure  5  shows  the  comparison  between  the  exper¬ 
imental  and  corresponding  predictions  using  Eq.  (  I )  The 
agreement  between  the  two  is  excellent  and  not  only  confirms 
the  validity  of  Eq .  ( 1 )  in  predicting  the  thermal  expansion  coef¬ 
ficients  of  these  composites,  but  also  indirectly  confirms  the 
value  of  6. 5  x  It)  ',<5C  1  measured  for  the  fibers. 

IV.  Summary  and  Conclusions 

( 1 )  The  thermal  expansion  coefficients  of  SiC  monofila¬ 
ments  were  measured  to  be  6.5  ±0.5  x  10' ""C  'in  the  tem¬ 
perature  range  200°-900°C. 

(2)  The  linear  theory  predicting  the  CTE  for  unidirectional 
fiber-reinforced  ceramic  matrix  composites  was  shown  to  be 
satisfactory  and  indirectly  confirm.-,  the  value  of  the  CTE  mea¬ 
sured  for  the  fibers. 

(3)  Heat  treatment  of  fhe  fibers  ( I400°C.  90  min)  resulted 
in  a  noticeable  increase  in  the  intensity  and  a  sharpening  of  the 
(3-SiC  peak,  indicating  some  grain  growth. 

(4)  The  reason  for  the  anomalously  high  CTE  measured  for 
the  fibers  is  not  entirely  clear,  but  the  XRD  results  seem  to  indi¬ 
cate  that,  possibly,  the  presence  of  an  amorphous  grain-bound¬ 
ary  phase  in  the  fibers  with  a  higher  CTE  than  bulk  SiC  could 
be  responsible. 
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